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ABSTRACT 


A  FORTRAN  Control  Program  for  Engineering  Synthesis  (COPES)  is 
presented  for  solving  engineering  design  problems.  This  program  uses 
the  general-purpose  optimization  code  ADS  (Automated  Design  Synthesis) 
to  provide  a  versatile  design  optimization  capability.  The  Pr°8ra® 
minimizes  or  maximizes  a  numerically  defined  objective  function  subject 
to  a  set  of  constraint  limits.  In  addition  to  direct  optimization, 
COPES  includes  approximation  techniques  as  well  as  the  capability  o 
perform  conventional  parametric  studies.  Programming  guidelines  are 
presented,  followed  by  sample  input  data  and  output  for  the  various 
program  options. 
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1.0  INTRODUCTION 


Most  design  processes  require  the  minimization  or  maximization  of 
some  parameter  which  may  be  called  the  design  objective.  For  the  design 
to  be  acceptable,  it  must  satisfy  a  variety  of  limitations  which  are 
referred  to  as  constraints.  While  much  of  the  practical  design  task 
cannot  be  automated,  the  engineering  tasks  which  are  quantifiable  are 
usually  well  suited  for  computer  solution. 

To  the  extend  that  the  problem  can  be  stated  in  rumerical  terms,  a 
computer  program  can  be  written  to  perform  the  necessary  calculations. 
For  this  reason,  computer  analysis  is  commonplace  in  most  engineering 
organizations.  For  example,  in  structural  design,  the  materials,  loads, 
and  configuration  may  be  defined,  and  a  finite  element  analysis  computer 
code  is  used  to  calculate  stresses,  deflections,  frequencies  and  other 
response  quantities  of  interest.  If  any  of  these  parameters  are  not 
within  the  prescribed  bounds,  the  engineer  may  change  the  structural 
member  sizes  or  geometry^  and  re-run  the  program.  The  computer  code 
therefore  provides  only  the  analysis  of  a  proposed  design,  with  the 
engineer  making  the  actual  design  decisions.  This  approach  to  design, 
which  may  be  called  computer-aided  design  is  commonly  used  today. 

Another  common  use  of  analysis  codes  is  in  parametric  studies.  For 
example,  an  aircraft  trafectory  analysis  code  may  be  run  repetitively 
for  several  payloads,  calculating  the  aircraft  range,  to  determine  the 
range-payload  dependency. 

A  logical  extension  to  computer-aided  design  is  fully  automated 
design,  where  the  computer  makes  the  actual  design  decisions,  or  to 
perform  parametric  studies  with  mlnimium  man-machine  interaction.  The 
purpose  of  the  COPES  program  is  to  provide  this  automated  design  and 
parametric  capability.  The  user  must  provide  a  FORTRAN  program,  in 
subroutine  form,  for  analysis  of  the  particular  problem  being  solved. 
This  analysis  program  is  written  according  to  a  simple  set  of  guidelines 
so  that  it  can  be  easily  coupled  to  the  COPES  program  for  automated 
design  synthesis. 

The  program  described  here  is  a  modification  of  the  earlier  version 
of  the  COPES  program  which  used  the  CONMIN  optimization  program  (refs. 
1,  2).  The  principal  differences  are  in  the  use  of  the  ADS  optimization 
program  (refs.  3,  4),  deletion  of  the  two-variable  function  space 
option,  optimum  sensitivity,  and  minor  changes  in  the  data  input. 

This  document  describes  the  capabilities  of  the  COPES  program  and 
its  usage.  A  simple  design  example  is  first  presented  to  demonstrate 
the  program  capabilities.  Guidelines  are  given  for  writing  analysis 
codes  which  can  be  coupled  directly  to  COPES.  Finally  the  data 
organization  is  outlined  and  sample  data  is  presented.  It  is  assumed 
that  the  user  of  COPES  is  familiar  with  the  AOS  user's  manual  (ref.  4), 
or  at  least  has  that  manual  available  for  ready  reference. 

The  most  commonly  used  option  of  COPES  is  for  design  optimization. 
Two  approaches  are  available  for  this  purpose.  The  first  (NCALC*2)  is 
direct  optimization  of  the  function  by  the  ADS  optimization  sub¬ 
program.  An  alternative  to  this  is  through  the  use  of  approximation 
techniques  (NCALCa6,  ref.  3).  This  second  option  is  often  more 
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efficient  for  problems  of  under  six  design  variables,  but  which  require 
costly  analysis,  especially  if  multiple  optimizations  are  performed. 
Also,  this  option  is  useful  for  design  where  the  analysis  Information  is 
not  actually  provided  by  a  FORTRAN  program,  but  is  instead  the  result  of 
experimentation  (ref.  6). 

In  the  following  section,  a  simple  example  using  the  classical 
three-bar  truss  (ref.  7)  is  given  to  demonstrate  problem  formulation  and 
solution  using  COPES /ADS. 


1.1  Design  Example 


Assume  it  is  required  to  design  the  simple  planar  truss  shown  in 
Figure  1.  The  objective  is  to  find  the  minimum  volume  of  material  which 
will  support  the  two  (independent)  loading  conditions  shown. 


Figure  l:  Three-Bar  Truss 


The  design  problem  is  stated  as  follows:  Find  the  cross-sectional 
areas  Ai  and  A2  (with  A3  -  Aj  for  symmetry)  which  will 


Minimize 

VOL  -  UIV2A!  +  A2  +  V2A3] 

(1) 

Subject  to; 

-15,000  <  CTjl j  ^  20,000  i-1,3 

j-1,2 

(2) 

At  >  0.01  i-1,3 

A3  -  Ai 


where  is  the  stress  in  member  i  under  load  condition  j. 


(3) 

(4) 


The  equations  necessary  to  calculate  the  member  stresses  are  as 
follows: 


Loading  Condition  1 


A2A3P1 


H  A^  A^ (A1A2+A2A3+  2AiA3; 
2A3pl 

cr  -  - - 

A1A2+A2A3+  2A]A3 

t., - - 

A1A2+A2A3+  2A]A3 


Loading  Condition  2_ 

„  “A2P1 

<Ti  2  *  . . 

AjA2+A2A3+  2A]A3 

<*22  “  _ _ 

A1A2+A2A3+  2A1A3 

a  p2  A1A2P2 

^  A^  A3  ( A2+A2A3+  2A1A3) 


Equation  (l)  is  the  objective  function  (volume)  and  Equations  (2) 
are  constraints  on  stress  in  the  members.  Equation  (3)  represents  side 
constraints"  on  the  members  to  prevent  the  optimizer  from  proposing  a 
design  with  negative  member  areas.  Equation  (4)  is  a  design  variable 
"linking"  equation  to  maintain  symmetry.  The  units  on  dimensions, 
stresses  and  loads  are  assumed  to  be  consistent  and  the  units  used  here 
are  understood  to  be  inches  and  pounds. 

This  design  problem  may  be  stated  in  standard  mathematical  form  for 
optimization  as: 

Minimize  F(X)  (5) 

Subject  to;  Gj(X)  <  0  j-I,m  (6) 

x{  <  Xi  <  Xi  i-l,n  (7) 

where  X  is  a  vector  containing  the  independent  design  variables,  A^  and 
A2  and  the  Gj(X)  are  the  constraints  defined  by  Equations  (2).  The 
objective  ana  constraints  are  functions  of  the  design  variables,  but 
need  not  be  explicit  as  given  here.  These  functions  can  be  evaluated  by 
any  means  at  our  disposal  including  implicit  iteration  techniques  or 
experimentation.  It  is  assumed,  however,  that  the  functions  are 
continuous  and  are  defined  everywhere  in  the  range  between  X*  and  X^« 
The  bounds  on  the  design  variables  X  given  by  Equation  (7)  are  referred 
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Co  as  side  constraints  and  limit  the  region  of  search  for  the  optimum. 
If  any  side  constraint  is  not  needed,  its* value  should  be  set  to  a  large 
negative  or  positive  value  for  Xj  or  x“,  respectively.  While  the  three- 
bar  truss  problem  given  above  is  straightforward,  its  solution  is  not 
trivial  because  it  is  not  known  which  constraints  will  be  critical  [i.e. 
which  Gj(X)  ■  0]  at  the  optimum. 

One  design  approach  is  to  assume  that  each  member  will  be  stressed 
to  its  limit  under  at  least  one  of  the  loading  conditions.  The  stresses 
are  calculated  and  each  member  is  resized  to  its  allowable  stress.  Thi3 
must  be  done  iteratively  because  the  stresses  are  functions  of  the 


no  part  in  the  design  decisions  and  so  there  is  no  assurance  that  this 
is  actually  a  minimum  volume  design.  More  importantly,  it  is  desirable 
to  devise  techniques  for  optimum  design  of  systems  which  may  be  defined 
by  more  than  two  variables  and  by  much  more  complex  analysis,  where  the 
optimum  cannot  easily  be  found  by  trial  and  error. 

The  COPES  program  provides  this  general  capability  by  the  use  of 
the  optimization  program  ADS.  To  use  this  design  capability,  a  FORTRAN 
code  must  be  provided  which  will  calculate  the  various  parameters.  In 
writing  the  analysis  code:  1)  it  is  written  in  subroutine  form  with 
SUBROUTINE  ANALIZ  (ICALC)  as  the  main  routine,  2)  it  is  segmented  into 
INPUT,  EXECUTION  and  OUTPUT  and  3)  all  parameters  which  may  be  design 
variables,  objective  functions  or  constraints  are  contained  in  a  single 
labeled  common  block  called  GLOBCM.  The  principal  motivation  for  the 
COPES  program  is  that  the  user  does  not  have  to  identify  the  design 
variables,  objective  or  constraints  in  advance  and  that  the  user  does 
not  have  to  formulate  the  problem  in  the  standard  form  for  optimization 
given  by  Equations  (5-7).  The  design  parameters  are  defined  by  data  at 
the  time  of  execution  and  the  analysis  information  is  converted  by  COPES 
to  the  standard  form  for  use  in  ADS.  Thus  the  COPES  program  greatly 
simplifies  the  use  of  optimization  techniques  for  routine  design  tasks. 

To  demonstrate  the  simplicity  with  which  a  design-oriented  analysis 
code  can  be  written,  the  following  FORTRAN  subrou  tine  was  written  for 
the  analysis  of  the  three-bar  truss  in  Figure  1.  Here  Pi  “  ?2  " 


'-T-1'.  J.‘  IJ.?  ..'A*  IJ.11  -.,p  ■  ■  .  f.  .r,  ■  . 

• 

.  *  .  *« 

SUBROUTINE  ANALIZ  (ICALC) 

•  i 

COMMON  /GLOB CM/  Al,A2,A3,VOL,SIGll,SIG21,SIG31,SIG12,SIG22,SIG32 

.  ■  _  “  * 

IF  ( ICALC.  GT.l)  GO  TO  10 

.*  ‘ 

C  -  INPUT  OR  INITIALIZATION. 

Al-1. 

• 

A2-1. 

*  V  V 

A3-1  • 

•  •, 

E-10.0E+6 

- 

H-10. 

a  '  • 

P-20000. 

a*  * 

RETURN 

• 

C  -  EXECUTION 

10  CONTINUE 

.*  "l 

VOL-H* ( (A1+A3 )*SQRT( 2 . )+A2 ) 

DEN0M-Al*A2+A2*A3+SQRT(2. )*A1*A3 

SIG11-P /A1-A2*A3*P / ( A1*DEN0M  ) 

- .  1 

SIG21-SQRT<2. )*A3*P/DENCW 

• 

f  .« 

SIG31— A2*P/DEN0M 

'  a 

SIG12-SIG31 

SIG22-SQRT(2. )*A1*P/DEN0M 

SIG32-P/A3-A1*A2*P/(A3*DEN0M) 

-  *  . 

IF  (ICALC.LT. 3)  RETURN 

C  -  OUTPUT. 

ft. 

-  ‘ 

WRITE(6,20)  P,H,E 

20  FORMAT  ( 1H1 , 4X , 22HTHREE-BAR  TRUSS  OUTPUT / /5X, 3HH  -.E12.5/5X, 

1  3 HP  -,£12.5/5X,3HE  -,E12.5) 

WRITE (6, 30)  A1,A2,A3 

30  FORMAT  (/5X.4HA1  -,Ei2.5/5X,4HA2  -,E12.5/5X,4HA3  -,E12.5) 
WRITE(6,40)  SIGH ,SIG2I,SIG3I ,SIG12,SIG22,SIG32 
40  FORMAT  (/5X.7HSIG11  «,E12.5/5X,7HSIG21  -,E12.5/5X,7HSIG31  -, 

I  E12.5/5X.7HSIG12  -,E12.5/5X,7HSIG22  -,E12.5/5X,7HSIG32  -,E12.5) 
WRITE (6, 50)  VOL 

50  FORMAT  (/5X.8HV0LUME  -,E12.5) 

RETURN 

END 


This  routine  may  be  executed  as  an  analysis  program  by  the 
following  main  program: 

C  MAIN  PROGRAM  TO  EXECUTE  SUBROUTINE  ANALIZ. 

DO  10  1-1,3, 2 
10  CALL  ANALIZ  (I) 

STOP 

END 

Moreover,  ANALIZ  can  be  coupled  directly  to  the  COPES  program  to 
perform  this  same  function,  or  to  perform  optimization  or  trade-off 
studies. 


---•I 
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This  subroutine  was  coupled  to  COPES  and  an  optimization  was 
performed  (NCALC  ■  2)  to  give  the  following  design: 

THREE -BAR  TRUSS 

P  -  20000 
H  -  10 
E  -  1.0E+7 

A1  -  0.7847 
A2  -  0.4201 
A3  -  0.7847 

SIGH  -  19996 
SIG21  -  14506 
SIG31  -  -5490 
SIG12  -  -5490 
SIG22  *  14506 
S1G32  -  19996 

VOLUME  -  26.396 

The  parameters,  P,  H  and  E  were  input  to  the  program  and  the 
remaining  parameters  were  calculated.  The  independent  design  variables, 
A1  and  A2  were  changed  during  the  optimization  and  symmetry  was  imposed 
so  that  A3  *  Al.  Note  that  A3  could  have  been  treated  as  an  independent 
variable  as  well  but,  assuming  symmetry  of  the  structure  was  desired, 
this  “linking”  of  the  design  variables  simplifies  the  optimization  task. 

This  design  was  achieved  with  21  calls  to  ANALIZ  with  ICALC  »  2 
(that  is,  the  design  process  required  21  separate  analyses).  This 
design  could  surely  have  been  found  with  fewer  analyses  were  it 
performed  by  hand  calculations.  However,  once  having  written  the 
analysis  subroutine,  numerous  other  designs  may  be  obtained  for 
different  materials,  loading,  or  design  stresses  with  only  minimal 
effort.  Furthermore,  the  design  obtained  here  with  a  volume  of  26.4 
in. 3  is  below  the  design  OF  38.2  in.3  found  by  ratloing  the  members 
according  to  their  stresses  as  may  be  done  in  traditional  design 
methods. 

This  simple  design  example  underscores  the  power  of  numerical 
optimization  techniques  and  the  ease  with  which  they  may  be  applied. 
The  key  to  efficient  use  ot  the  COPES  program  is  the  requirement  that 
the  ANALIZ  code  be  written  in  a  standard  format  so  that  it  may  be 
coupled  to  COPES  without  modification.  The  following  section  contains 
guidelines  for  writing  a  design-oriented  analysis  code. 


2.0  PROCRAMMING  GUIDELINES 

la  developing  any  computer  code  for  engineering  analysis,  it  is 
prudent  to  write  the  code  in  such  a  way  that  it  is  easily  coupled  to  a 
general  synthesis  program  such  as  COPES.  Therefore,  a  general 
programming  practice  is  outlined  here  which  in  no  way  inhibits  the  use 
of  the  computer  program  in  its  traditional  role  as  an  analytic  tool,  but 
allows  for  simple  adaptation  to  COPES.  This  approach  is  considered  good 
programming  practice  and  provides  considerable  flexibility  of  design 
options.  Only  five  basic  rules  oust  be  followr-*: 

1.  Write  the  code  in  subroutine  form  with  the  primary  routine  called 
as;  SUBROUTINE  ANALIZ  (ICALC).  The  name  ANALIZ  is  compatible  with 
the  COPES  program  and  ICALC  is  a  calculation  control.  Note  that 
subroutine  ANALIZ  may  call  numerous  other  subroutines  as  required  to 
perform  the  necessary  calculations. 

2.  Segment  the  program  into  INPUT,  EXECUTION  and  OUTPUT.  The 
calculation  control,  ICALC,  will  determine  the  portion  of  the 
analysis  code  to  be  executed. 

ICALC  -  1:  The  program  reads  all  data  required  to  perform  the 
analysis.  Also,  any  initialization  of  constants  which  will  be  used 
repetitively  during  execution  is  done  here.  This  initial  input 
information  is  printed  here  for  later  reference  and  for  program 
debugging . 

ICALC  »  2:  The  program  performs  the  execution  phase  of  the  analysis 
task.  No  data  reading  or  printing  is  done  here,  except  on  user** 
defined  intermediate  files.  Data  may  be  printed  here  during  program 
debugging,  in  which  case  it  should  be  controlled  by  a  print  control 
parameter  which  is  read  during  input.  In  this  way,  this  print  may 
be  turned  off  after  the  program  is  debugged,  but  may  be  used  again 
during  future  program  expansion  and  debugging.  THe  reason  that 
printing  is  not  normally  allowed  during  execution  is  that  when 
optimization  is  being  done,  the  code  will  be  called  many  times  with 
ICALC  ■  2,  resulting  in  voluminous  print. 

ICALC  -  3:  The  results  of  the  analysis  are  printed.  Also,  the 
essential  input  parameters  which  may  have  been  changed  during 
optimization  should  be  printed  here  for  easy  reference.  When 
ICALC  *  3,  it  is  good  practice  to  execute  the  analysis  also  because 
the  analysis  information  may  not  be  up  to  date.  This  is  because  the 
optimization  program  may  not  choose  the  last  design  considered  as 
the  optimum. 

In  Summary,  When: 

ICALC  -  1,  Read  input  data. 

ICALC  ■  2,  Execute  the  analysis. 

ICALC  ■  3,  Execute  the  analysis  and  print  the  results. 


3.  Store  all  parameters  which  may  be  design  variables,  objective 
functions  or  constraints  in  a  single  labeled  COMMON  block  called 
GLOBCM.  The  order  in  which  they  are  stored  is  arbitrary.  Also,  it 
is  not  necessary  to  identify  the  design  variables,  objective 
function  and  constraints  in  advance.  Any  parameters  which  may  be 
used  for  optimization  or  parametric  studies  may  be  stored  in  GLOBCM. 
A  listing  of  the  COPES  program  should  be  checked  to  see  how  many 
parameters  may  be  stored  in  GLOBCM  (this  is  the  dimensioned  size  of 
the  array  called  ARRAY).  Initial  distribution  of  COPES/ADS  allows 
for  500  parameters.  If  more  storage  is  needed,  this  array  can  be 
dimensioned  accordingly  in  COPES  and  the  parameter  NARRAY  at  the 
beginning  of  the  program  changed  to  the  new  dimensioned  size  of 
ARRAY. 

A.  During  execution  or  output,  no  parameters  which  are  read  during 

input  should  be  updated.  For  example,  if  variable  Q  is  initialized 
during  input,  the  execution  segment  must  not  update  Q  such  as 
Q-Q+3.2.  Instead,  a  new  variable,  Y*Q+3.2  should  be  defined. 

5.  Write  all  programs  in  standard  language,  avoiding  machine  dependent 
capabilities.  While  this  guideline  is  not  essential  to  the  use  of 
the  analysis  code  with  the  COPES  program,  it  makes  the  analysis  code 
much  more  transportable  between  different  computer  systems,  a 
capability  which  usually  justifies  a  slight  reduction  in  efficiency 
on  a  given  machine. 

Adherence  to  these  guidelines  not  only  leads  to  a  more  readable  and 
machine  independent  computer  code,  but  allows  this  code  to  be  coupled  to 
the  COPES  program  without  modification. 

Having  written  the  analysis  code,  it  may  be  executed  either  with  a 
simple  main  program  or  within  the  COPES  program  to  perform  analysis.  To 
insure  that  guideline  A  is  followed,  the  main  test  program  given  belcw 
is  recommended.  Note  that  this  program  calls  ANALIZ  twice  with  ICALC“3 
(assuming  execution  is  also  performed  when  ICALC*3)  to  show  that  the 
same  result  is  obtained  repetitively. 

C  MAIN  PROGRAM  TO  CHECK  SUBROUTINE  ANALIZ. 

C  READ,  EXECUTE  AND  PRINT. 

DO  10  ICALC-1,3,2 
10  CALL  ANALIZ  (ICALC) 

C  EXECUTE  AND  PRINT  AGAIN  TO  BE  SURE  THE  RESULTS  ARE  THE  SAME. 

DO  20  ICALC-2,3,2 
20  CALL  ANALIZ  (ICALC) 

STOP 

END 

Once  the  analysis  code  has  been  written,  it  can  be  coupled  to  the 
COPES  program  without  modification.  If  it  is  desired  to  perform  a 
simple  analysis  using  COPES,  only  three  data  cards  are  required  for  the 
COPES  program,  namely  the  TITLE  card,  the  control  parameter,  NCALC  -  l, 
and  an  END  card.  If  the  optimization  or  parametric  analysis 
capabilities  of  COPES  are  to  be  used,  additional  data  must  be  read. 
This  data  will  identify  which  parameters  in  the  global  common  block, 
GLOBCM,  are  used. 


3.0  COPES  DATA  MANAGEMENT 


la  order  to  perform  design  operations,  the  COPES  program  must 
access  the  data  in  common  block  GLOBCM.  This  is  done  by  defining  the 
location  in  GLOBCM  where  a  specified  parameter  resides.  For  example,  * _ 

consider  the  common  block  for  the  three-bar  truss:  i 

COMMON  /GLOBCM/  A1,A2,A3, VOL, SIG11, SIG21,SIG31,SIG12,SIG22,SIG32 

The  volume  of  material,  VOL,  is  the  fourth  parameter  in  the  common 
block;  that  is,  it  resides  in  location  4,  referred  to  as  the  global 
location  number.  Similarly  the  stress  in  bar  1  under  loading  condition  9 

2,  SIG12,  is  in  global  location  8.  Thus,  the  parameters  are  referred  to 
by  their  respective  location  numbers  in  global  common. 

For  convenience  in  preparing  data  for  the  COPES  program,  a  simple 
"CATALOG**  of  parameters  may  be  defined.  For  the  three-bar  truss,  this 
catalog  would  be:  - 


GLOBAL 

LOCATION 

FORTRAN 

NAME 

DEFINITION 

1 

A1 

Cross-sectional  area  of  member  1. 

2 

A2 

Cross-sectional  area  of  member  2. 

3 

A3 

Cross-sectional  area  of  member  3. 

4 

VOL 

Volume  of  material  in  the  structure. 

5 

SIGH 

Stress  in  member  1  under  load  condition 

1. 

6 

S1G21 

Stress  in  member  2  under  load  condition 

1. 

7 

SIG31 

Stress  in  member  3  under  load  condition 

1. 

8 

SIG12 

Stress  in  member  1  under  load  condition 

2. 

9 

SIG22 

Stress  in  member  2  under  load  condition 

2. 

10 

SIG32 

Stress  in  member  3  under  load  condition 

2. 

As  another  example,  consider  a  global  common  block  containing 
arrays : 

COMMON  /GLOBCM/  A,Y(10) ,Q,C(2,2) ,H 
The  variable  catalog  for  this  common  block  is: 


GLOBAL 

LOCATION 

FORTRAN 

NAME 

DEFINITION 

1 

A 

Area 

2 

*<10) 

Vector 

of  Y-coordinates 

12 

Q 

e 

13 

C(2,2) 

• 

17 

H 

etc. 

The  dimensions  are  given  with  the  FORTRAN  name  as  a  reminder  that 
the  parameter  is  an  array.  In  this  case,  the  third  parameter  in  the  Y- 
array  is  in  global  location  4.  Remembering  that  arrays  are  stored 


column  by  column,  che  C(l,2)  array  location  is  in  global  location  IS. 

It  will  be  seen  that  identif ying-parameters  according  to  their 
location  in  GLOBCM  provides  a  great  deal  of  flexibility  in  using  the 
COPES  program  for  design. 


4.0  COPES  TERMINOLOGY 

The  COPES  program  provides  four  specific  capabilities: 

1.  Simple  analysis,  just  as  if  COPES  was  not  used. 

2.  Optimization  -  Minimization  or  maximization  of  one  calculated 
parameter  with  limits  imposed  on  other  parameters. 

3.  Approximation  techniques  -  Optimization  using  approximations.  Often 
more  efficient  than  standard  optimization  for  up  to  six  design 
variables  or  if  multiple  optimizations  are  to  be  performed. 

4.  Parametric  analysis  -  The  effect  of  changing  one  or  more  design 
variables  on  one  or  more  calculated  parameters. 

In  defining  the  data  required  to  execute  the  COPES  program,  the 
following  definitions  are  useful: 

Design  Variables  -  Those  parameters  which  the  optimization  program  is 
allowed  to  change  in  order  to  improve  the  design.  Design  variables 
appear  only  on  the  right  hand  side  of  equations  in  the  analysis  program. 
COPES  considers  two  types  of  design  variables,  independent  and 
dependent.  If  two  or  more  variables  are  always  required  to  have  the 
same  value  or  be  in  a  constant  ratio,  one  is  the  independent  variable 
while  the  remaining  are  dependent  variables.  For  example,  in  the  three- 
bar  truss,  A1  may  be  the  independent  variable,  while  A3-A1  is  a 
dependent  variable,  treated  this  way  to  insure  symmetry  in  the  design. 

Objective  Function  -  The  parameter  which  is  to  be  minimized  or  maximized 
during  optimization.  Objective  functions  always  appear  on  the  left  side 
of  equations,  unless  the  objective  function  is  also  a  design  variable 
(i.e.  the  cross-sectional  area  of  member  A1  of  the  three  bar  truss  may 
be  minimized  as  the  objective  if  there  is  a  particular  reason  to  do  so. 
In  that  case,  A1  would  be  both  a  design  variable  and  the  objective 
function).  An  objective  function  may  be  linear  or  nonlinear,  implicit 
or  explicit,  but  must  be  a  function  of  the  design  variables  to  be 
meaningful. 

Constraint  -  Any  parameter  which  must  not  exceed  specified  bounds  for 
the  design  to  be  acceptable.  Constraint  functions  always  appear  on  the 
left  side  of  equations.  Just  as  for  objective  functions,  constraints 
may  be  linear  or  nonlinear,  implicit  or  explicit,  but  must  be  functions 
of  the  design  variables. 

Constraint  Set  -  A  group  of  constraints  which  appear  consecutively  in 
the  global  common  block,  GLOBCM,  and  which  all  have  the  same  limits 
imposed.  This  is  a  convenience  which  allows  several  constraints  to  be 
identified  with  a  minimum  of  data. 


Global  Comaon  -  Common  block  GLOBCM  containing  design  Information* 

Global  Location  -  Location  of  a  particular  parameter  In  GLOBCM. 

Strategy  -  Optimization  strategy  such  as  Sequential  Quadratic 
Programming  . 

Optimizer  -  Actual  algorithm  to  perform  the  optimization. 

One~Dlmenaional  Search  ~  Line  search  routine  used  by  the  Optimizer. 

5.0  COPES  DATA 

The  COPES  program  reads  data  from  unit  5  and  writes  output  to  unit 
6.  Units  10  and  11  are  used  as  scratch  files.  The  scratch  file  numbers 
may  be  changed  by  changing  two  cards  at  the  beginning  of  the  COPES 
program.  The  FORTRAN  names  defining  the  file  numbers  are  ISCR1  and 
ISCR2. 

In  order  to  execute  the  COPES  program,  it  is  necessary  to  provide 
data  for  COPES,  followed  by  data  for  the  ANALIZ  program  which  is  coupled 
to  COPES.  This  section  defines  the  data  which  is  required  by  COPES. 
The  data  is  segmented  into  "BLOCKS"  for  convenience.  All  formats  are 
alpha-numeric  for  TITLE  and  END  cards,  F10  for  real  data  and  110  for 
integer  data.  Unformated  input  is  allowed  for  COPES  data. 

The  COPES  data  begins  with  a  TITLE  card  and  ends  with  an  END  card. 
This  is  followed  by  data  to  be  read  by  the  user  supplied  subroutine 
ANALIZ  when  ICALC  -  1. 

Comment  cards  may  be  inserted  anywhere  in  the  COPES  data  stack 
prior  to  the  END  card,  and  are  identified  by  a  dollar  sign  ($)  in  column 
1  of  the  card. 

5.1'  Unformatted  Data  Input 

While  the  user's  sheet  defines  COPES  data  in  formatted  fields  of 
ten,  the  data  may  actually  be  read  in  a  simplified  fashion  by  separating 
data  by  commas  or  by  one  or  more  blanks.  If  more  than  one  number  is 
contained  on  an  unformatted  data  card,  a  comma  must  appear  somewhere  on 
the  card.  If  exponential  numbers  such  as  2.5+10  or  2.5E+10  (depending 
on  the  computer  system)  are  read  on  an  unformatted  card,  there  must  be 
no  embedded  blanks.  Unformatted  cards  may  be  intermingled  with 
formatted  cards.  Real  numbers  on  an  unformatted  card  should  have  a 
decimal  point. 

Examples 

Unformatted  data: 

5, 7, 1.3, 1.0+20,0, -5.1 
5, 7, 1.3,1. 0+-20  ,,-5.1 
5  7  1.3  1.0+20,,  -5.1 
5  7  1.3,  1.G+-20  0  -5.1 


5 


Unformatted  data 


2 

2,3 
2  3 

i 

Equivalent  formatted  data: 

Col  10  20  30 

2  Note.  This  data  has  been  right  justified. 

2  3 

2  3  Note.  This  data  contains  no  commas,  so  it 

is  assumed  to  be  formatted  already. 


Unformatted  data: 
1,2,3,4,5,6,7,8,9,10,11 
Equivalent  formatted  data: 


5.2  Optimization  Example 

For  the  three-bar  truss  example,  it  is  now  a  simple  matter  to 
create  the  COPES  data  to  solve  the  design  problem.  An  essential  part  of 
this  is  problem  formulation.  Therefore,  the  basic  design  questions  are 
first  presented,  followed  by  the  necessary  data  to  solve  this  problem 
using  COPES /ADS. 

What  Must  We  Knew? 

Number  of  Independent  design  variables: 

Al,  A2  (A3-A1)  NDV-2 

4 

*—  Dependent  Variable 
Total  number  of  design  variables: 

Al,  A2,  A3  NDVTOT-3 

What  are  the  locations  of  the  design  variables? 

1,  2  and  3 

What  is  the  objective? 

Minimize  VOL 

^-Location  4 

SGNOPT*— 1 . 0  (+1.0  would  maximize) 

Constraints? 

Yes,  -15,000  <  SIGU  <  20,000 


1  Set,  Locations  5-10,  Nonlinear  inequality 


COPES  Data  to  Perform  Optimization 


The  data  is  simply  listed  here.  Further  definitions  of  the 
parameters  are  given  in  the  data  block  input  in  the  following  section. 

$  BLOCK  A:  TITLE 

3-BAR  TRUSS  DESIGN  EXAMPLE 

$  BLOCK  B: 

2.2  Optimization,  2  independent  variables 
$  BLOCK  C: 

0,4,8,2020  Strategy,  Optimizer,  One-Dimensional  Search, 

Print  control 

$  BLOCK  D:  OMIT 
S  BLOCK  E:  CHIT 
$  BLOCK  F: 

3, 4, -1.0  NDVTOT ,  IOBJ,  SGNOPT 

S  BLOCK  G: 

0. l,1.0E+20  Design  Variable  Bounds 

0.1,i.0E+20 

$  BLOCK  H: 

1.1  Design  Variable  Number  and  Global  Location 

2.2 

1.3 

S  BLOCK  I: 

l  Number  of  Constraint  Sets 

$  BLOCK  J: 

5,10  Global  Locations  of  Constrained  Variables 

-15000. ,0. ,20000.  Lower  and  Upper  Bounds 

$  BLOCK  R: 

END  Required  END  card 

This  data  was  used,  together  with  the  ANALIZ  subroutine  in  Section 
1.1  to  create  the  design  given  there.  This  example  is  used  in  Section 
6.0  to  demonstrate  all  basic  options  of  COPES  and  sample  output  is  given 


5.3  Copes  Data  Input 

The  necessary  input  data  to  the  COPES  program  is  given  here.  Note 
that,  depending  on  the  program  options  to  be  used,  much  of  this  data  may 
be  omitted. 


DATA  BLOCK  A 

DESCRIPTION;  Title  Card. 

FORMAT  AND  EXAMPLE 

1 _ 2 _ 3 _ 4 _ 5 _ 6 _ 7 

TITLE 

THREE-BAR  TRUSS 

FIELD  CONTENTS 


8  FORMAT 


1-8 


Any  title  may  be  given  on  this  card,  up  to  80  characters 


DATA  BLOCK  B 


DESCRIPTION:  Program  Control  Parameters. 
FORMAT  AND  EXAMPLE 


123456  FORMAT 


NCALC 

NDV 

NAPRX 

NPV 

IPNPUT 

IPDBG 

6110 

2 

2 

1 

3 

0 

0 

FIELD 

1  NCALC: 

0  - 

1  - 

2  - 

3  - 

4  - 

2  NDV : 

3  NAPRX: 

0  - 
1  - 

4  NPV: 

5  IPNPUT : 

0  - 

1  - 
2  - 

6  IPDBG: 
REMARKS 


CONTENTS 


Calculation  Control 

Read  input  and  stop.  Data  of  blocks  A,  B  and  R  Is 
required.  Remaining  data  is  optional. 

One  cycle  through  the  program.  The  same  as  executing 
ANALIZ  as  a  stand-alone  program.  Data  of  blocks  A,  B 
and  R  is  required.  Remaining  data  is  optional. 
Optimization.  Data  of  blocks  A-J  and  R  is  required. 
Remaining  data  is  optional. 

Optimization  Using  Approximation  Techniques.  Data  of 
blocks  A-0  and  V  is  required.  Remaining  data  is 
optional.  . 

Parametric  Study.  Data  of  blocks  A,  B  and  P-V  is 
required.  Remaining  data  is  optional. 

Number  of  independent  design  variables  for  optimization. 

Parameter  to  indicate  if  data  is  to  be  read  for 
optimization  using  approximation  techniques. 

No  data  will  be  read. 

Data  will  be  read. 

Number  of  parameters  to  be  varied  in  a  parametric  study. 
Input  data  print  control. 

Print  card  images  of  data  plus  formatted  print  of  input 
information. 

Formatted  print  of  input  information  only. 

No  print  of  input  data. 

Debug  print  control. 


1)  Data  to  be  read  in  the  following  blocks  depends  on  the  values  of 
NDV,  NAPRX  and  NPV.  Actual  program  execution  is  controlled  by  the  value 
of  NCALC,  and  it  is  assumed  by  COPES  that  the  necessary  data  has  been 
provided. 


DATA  BLOCK  C:  OMIT  IF  NDV  -  0  IN  BLOCK  B 

DESCRIPTION :  Integer  optimization  control  parameters. 

FORMAT  AND  EXAMPLE 


1 

2 

3 

4 

5 

6 

FORMAT 

ISTRAT 

10  PT 

IONED 

IPRINT 

IOVRDR 

IOVRDI 

6110 

0 

5 

7 

2020 

0 

0 

FIELD  CONTENTS 

1  Strategy  option  to  be  used.  DEFAULT  ■  0. 

2  Optimizer  option  to  be  used.  DEFAULT  ■  5. 

3  One-Dimensional  search  option  to  be  used.  DEFAULT  *  7. 

4  ADS  Print  control. 

5  Number  of  real  ADS  parameters  to  be  over-ridden.  Usually 
IOVRDR  ■  0  is  adequate  so  the  ADS  default  parameters  will  be 
used. 

6  Number  of  integer  ADS  parameters  to  be  over-ridden.  Usually 
IOVRDI  »  0  is  adequate  so  the  ADS  default  parameters  will  be 
used. 


REMARKS 

1)  See  Appendix  A  for  available  options  for  ISTRAT,  I0PT  and  IONED, 
as  well  as  guidelines  for  choosing  these  parameters. 

2)  See  Appendix  B  for  ADS  internal  parameters  which  may  be  modified 

by  the  user  to  "fine  tune"  the  optimization  to  a  specific  class 
of  problems. 


3)  See  the  ADS  User's  Manual  (ref.  4)  for  further  details. 


DATA  BLOCK  D:  OMIT  IF  IQVRDR  -  Q  IN  BLOCK  C 

DESCRIPTION ;  Identifiers  and  values  of  real  parameters  to  be  input. 
FORMAT  AND  EXAMPLE 


1  2  FORMAT 


ILOC 

VALUE 

110, F10 

11 

0.002 

FIELD  CONTENTS 

1  ILOC:  Location  in  array  WK  where  this  parameter  is  stored. 

2  VALUE:  New  value  of  the  parameter. 


REMARKS 

1)  See  Appendix  B  for  parameter  default  values,  definitions  and 
locations. 


i.i 


DATA  BLOCK  F :  OMIT  IF  NDV  -  0  IN  BLOCK  B 


DESCRIPTION:  Total  lumber  of  design  variables,  design  objective  global 
location  and  sign* 

FORMAT  AND  EXAMPLE 


1 

2 

3 

FORMAT 

NDV TOT 

IOBJ 

SGNOPT 

2I10,F10 

3 

4 

■a 

FIELD  CONTENTS 

1  NDV TOT:  Total  number  of  design  variables,  including  both 

dependent  and  independent  variables*  This  option  allows 
two  or  more  parameters  to  be  assigned  to  a  single  design 
variable*  The  value  of  each  parameter  is  the  value  of 
the  design  variable  times  a  multiplier,  which  may  be 
different  for  each  parameter*  DEFAULT  *  NDV. 

2  IOBJ:  Global  location  of  the  objective  function. 

3  SGNOPT:  Sign  used  to  identify  whether  the  objective  function  is 

to  be  maximized  or  minimized.  +1.0  indicates 
maximization.  -1*0  indicates  minimization.  If  SGNOPT 
is  not  unity  in  magnitude,  it  acts  as  a  multiplier  as 
well,  to  scale  the  magnitude  of  the  objective. 


REMARKS 

1)  If  maximization  is  to  be  performed,  the  objective  will  be 

changed  internally  and  the  objective  printed  by  ADS  will  appear 
opposite  in  sign.  This  is  because  ADS  minimizes  the  objective, 
so  if  maximization  is  performed,  ADS  will  minimize  the  negative 
of  the  true  objective. 


DATA  BLOCK  Gj_  OMIT  IF  NDV  -  0  IN  BLOCK  B 

DESCRIPTION :  Design  variable  bounds  and  initial  values. 

FORMAT  AND  EXAMPLE 


1 

2 

3 

FORMAT 

VLBI 

VUBI 

XI 

3F10 

0.01 

l.CH-20 

Bl 

NOTE:  READ  ONE  CARD  FOR  EACH  OF  THE  NDV  INDEPENDENT  DESIGN  VARIABLES. 
FIELD  CONTENTS 

1  VLBI:  Lover  bound  on  the  design  variable.  If  no  lower  bound  is 

to  be  imposed,  set  VLBI  to  a  large  negative  value. 

2  VUBI:  Upper  bound  on  the  design  variable.  If  no  upper  bound  is 

to  be  imposed,  set  VUBI  to  a  large  positive  value. 

3  XI:  Initial  value  of  the  design  variable.  Is  XI  is  non-zero, 

this  will  supercede  the  value  initialized  by  the  user- 
supplied  subroutine  ANALIZ. 
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DATA  BLOCK  H:  OMIT  IF  NDV  -  0  IN  BLOCK  B 
DESCRIPTION :  Design  variable  identification. 
FORMAT  AND  EXAMPLE 


2  3  FORMAT 


NOTE:  READ  ONE  CARD  FOR  EACH  OF  THE  NDVTOT  DESIGN  VARIABLES. 

FIELD  CONTENTS 

1  NDSGN:  Design  variable  number  associated  with  this  parameter. 

2  IDSGN:  Global  location  nuaber  associated  with  this  parameter. 

3  AMULT:  Constant  nultiplier  on  this  parameter.  The  value  of  the 

paremeter  will  be  the  value  of  the  design  parameter,  NDSGN, 
times  AMULT.  DEFAULT  -  1.0 

REMARKS 

1)  The  first  parameter  assigned  with  a  given  design  variable  is 
the  independent  variable.  Following  this,  any  parameters 
assigned  to  this  design  variable  are  dependent  variables. 

2)  The  value  of  AMULT  associated  with  an  independent  variable  should 
always  be  1.0. 


DATA  BLOCK  J:  OMIT  IF  SDV  ■  0  IN  BLOCK  B  OR  NCONS  -  0  IN  BLOCK  I 

DESCRIPTION:  Constraint  identification,  bounds  and  normalization 
factors. 

FORMAT  AND  EXAMPLE 


1 

2 

3 

4 

FORMAT 

ICON 

JCON 

LCON 

3110 

4 

8 

0 

BL 

FACT1 

BU 

FACT2 

4F10 

gggj 

wm 

20000. 

Da 

NOTE:  READ  TWO  CARDS  FOR  EACH  OF  THE  NCONS  CONSTRAINT  SETS. 

FIELD  CONTENTS 

1  ICON:  First  global  location  number  corresponding  to  this 

constraint  set. 

2  JCON:  Last  global  location  number  corresponding  to  this 

constraint  set.  DEFAULT  -  ICON. 


3  LCON:  Constraint  type  identifier. 

0,1  -  Nonlinear  Inequality  constraint. 

2  -  Linear  inequality  constraint. 

-1  -  Nonlinear  equality  constraint. 

-2  -  Linear  equality  constraint. 

FIELD  CONTENTS 

1  3L:  Lower  bound  on  the  constrained  variables. 

If  BL.LT.-1.0E+15,  the  lower  bound  is  ignored  by  COPES. 

2  FACT1:  Normalization  factor  on  the  lower  bound. 

DEFAULT  -  MAX  of  ABS(BL)  and  1.0. 

3  BU:  Upper  bajnd  on  the  constrained  variables. 

If  BU.GT.1.0E+15,  the  upper  bound  is  ignored  by  COPES. 


4 


FACT2:  Normalization  factor  on  the  upper  bound 
DEFAULT  -  MAX  of  ABS(BU)  and  1.0. 


REMARKS 


1)  The  normalization  factor  should  usually  be  defaulted. 

2)  The  constraint  functions  sent  to  AOS  are  of  the  fora: 

(BL  -  VALUE) /FACT 1  -LE.  0  and  (VALUE  -  BU)/FACT2  .LE.  0. 

3)  Each  constrained  parameter  is  converted  to  two  constraints  for 
AOS  unless  ABS(BL)  or  ABS(BU)  exceeds  1.QE+L5,  in  which  case  no 
constraint  is  created  for  that  bound. 

4)  In  ADS  Version  1.00,  equality  constraints  are  only  allowed  for 
ISTRAT  -  1-5. 
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DATA  BLOCK  Ki  OMIT  IP  NAPRX  -  0  IN  BLOCK  B 

DESCRIPTION :  Approximate  analysis/optimization  control  paraneters. 
FORMAT  AND  EXAMPLE 


1 

2 

3 

4 

5 

6 

7 

8 

FORMAT 

IPAPRX 

NXI 

NXFI 

NXA 

INCM 

ISCRX 

ISCRXF 

7110 

1 

4 

0 

1 

0 

0 

0 

KMIN 

KMAX 

NPMAX 

JNCM 

|jp|M 

5110 

0 

0 

0 

C 

0 

FIELD 


CONTENTS 


1  IPAPRX:  Print  control*  Values  of  0-4  with  Increasing  amounts  of 

print  for  larger  IPAPRX. 

2  NXI:  Number  of  candidate  X-vectors  read  as  data. 

3  NXFI:  Number  of  X-F  pairs  read  as  data. 

4  NXA:  0  or  1.  If  NXA  ■  1,  the  design  variables  read  by 

Subroutine  ANALIZ  form  a  candidate  X-vector. 


5  INOM:  Nominal  X-vector.  DEFAULT  ■  Best  available. 

6  ISCRX:  Unit  from  which  the  NXI  X-vectors  are  read.  DEFAULT  ■  5. 

If  ISCRX  is  not  S  (the  default),  the  data  is  assumed  to 
be  in  binary  format. 

7  ISCRXF:  Unit  from  which  the  NXFI  X-F  pairs  are  read. 

DEFAULT  -  5.  If  ISCRXF  is  not  5  (the  default),  the  data 
is  assumed  to  be  in  binary  format. 


FIELD  CONTENTS 

1  KMIN:  Minimum  number  of  approximation  iterations. 

DEFAULT  »  2*NDV  +  2  -  NXI  -  NXFI  -  NXA. 

2  KMAX:  Maximum  lumber  of  approximation  iterations. 

DEFAULT  -  NDV  +  NDV(NDV+l)/2  +  5.  If  KMAX  -  -1,  only 
one  approximate  optimization  is  performed  and  the 
resulting  proposed  optimum  is  printed,  together  with  the 
approximating  functions.  This  is  useful  when  the 
precise  functions  are  evaluated  experimentally  or  by 
another  program,  and  read  as  input  to  COPES  for  the  next 
approximate  optimization. 
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NPMAX:  Maximum  number  of  designs  retained  for  approximations. 
DEFAULT  -  ALL. 


0 

ir-,% 


4  JNOM:  Number  of  designs  after  which  the  best  design  is  picked 

as  nominal.  Before  this  the  most  recent  design  is 
nominal.  DEFAULT  -  2*NDV  +  2 

5  MAXTBM:  Terms  retained  in  the  Taylor  series  approximation. 

1  -  Linear  terms  only. 

2  -  Linear  terms  plus  diagonal  elements  of  the  Hessian 

matrix. 

3  -  Full  second-order  expansion. 

DEFAULT  -  3. 


► 


DATA  BLOCK  L:  OMIT  IF  NAPRX  -  0  IN  BLOCK  B 
DESCRIPTION :  Sounds  for  approximate  optimization. 

FORMAT  AND  EXAMPLE 

2345678  FORMAT 


DATA  BLOCK  M:  OMIT  If  NAPRX  -  0  IN  BLOCK  B 

DESCRIPTION :  Multipliers  on  bounds  for  approximate  optimization. 

FORMAT  AND  EXAMPLE 

1  2  FORMAT 


FIELD  CONTENTS 

1  XFACT1:  Multiplier  on  DXI  when  the  diagonal  elements  of  the 
Hessian  matrix  are  available.  DEFAULT  *  1.5. 


2  XFACT2:  Multiplier  on  DXI  when  all  elements  of  the  Hessian 
matrix  are  available.  DEFAULT  ■  2.0. 

REMARKS 

1)  These  factors  increase  the  move  limits  given  in  Block  L,  based  on 
the  assumption  that  larger  move  limits  are  possible  when 
nonlinear  information  is  available. 


DATA  BLOCK  N:  OMIT  IF  NAPRX  -  0  IN  BLOCK  B  OR  NXI  ■  0  IN  BLOCK  K 
DESCRIPTION :  Candidate  X-Vectors  for  approximate  optimization. 
FORMAT  AND  EXAMPLE 


12345678  FORMAT 


XII 

XI 2 

XI 4 

e  e  • 

•  e  e 

•  •  e 

•  •  • 

8F10 

0.8 

mu 

NOTE:  NXI  SETS  OF  DATA  ARE  READ  HERE. 

NOTE:  MORE  THAN  ONE  CARD  MAY  BE  READ  FOR  EACH  SET  OF  DATA. 
FIELD  CONTENTS 

1-8  XIJ:  Jth  value  if  the  Ith  X-vector,  J*1,NDV. 
REMARKS 


1)  If  ISCRX  is  not  5  (the  default  value),  data  read  in  this  block  is 
assumed  to  be  in  binary  format. 


DATA  BLOCK  OMIT  IF  NAPRX  -  £  IN  BLOCK  B  OR  NXFI  -  0  IN  BLOCK  K 
DESCRIPTION;  X-F  pairs  of  information  for  approximate  optimization. 
FORMAT  AND  EXAMPLE 


12345678  FORMAT 


XI 

X2 

X3 

X4 

•  e  e 

•  e  • 

•  •  • 

e  e  e 

8F10 

0.8 

0.6 

FI 

F2 

F3 

F4 

F5 

•  •  e 

e  e  e 

see 

8F10 

26.529 

20000. 

13474. 

-6525.7 

-6525.: 

NOTE:  NXFI  SETS  OF  DATA  ARE  READ  HERE. 

NOTE;  MORE  THAN  ONE  CAR D  MAY  BE  REQUIRED  FOR  XI  OR  FI. 

NOTE:  NDV  VALUES  OF  X  AND  NF  VALUES  OF  F  ARE  READ  FOR  EACH  SET  OF  DATA, 
WHERE  NF  -  1  +  THE  TOTAL  NUMBER  OF  CONSTRAINED  PARAMETERS 
(NOT  CONSTRAINT  SETS). 

FIELD  CONTENTS 

1  XI:  Ith  value  of  the  X-vector,  1*1, NDV. 

2  FI:  Ith  function  value,  I»1,NF. 

REMARKS 


1)  If  ISCRXF  is  not  5  (the  default  value),  data  read  in  this  block 
is  assumed  to  be  in  binary  format. 


DATA  BLOCK  P:  OMIT  IF  NPV  -  0  IN  BLOCK  B 
DESCRIPTION :  Parametric  variables* 

FORMAT  AND  EXAMPLE 

1  2  3  4  5 


NPF  IPRNTP 


8  FORMAT 


2110 


NPNI 

NPN2 

NPN3 

NPN4 

NOTE;  TWO  OR  MORE  CARDS  ARE  READ  HERE. 

FIELD  CONTENTS 

1  NPF:  Number  of  separate  functions  to  be  calculated  in  the 

parametric  study. 

2  IPRNTP:  Debug  print  control.  If  greater  than  zero,  detailed 

print  will  be  called  ats  each  step  in  the  parametric 
study.  DEFAULT  -  no  print. 

FIELD  CONTENTS 

1-8  NPNI:  Global  location  associated  with  the  I-th  function  to  be 
evaluated  in  the  parametric  study. 

REMARKS 

1)  More  than  eight  functions  are  allowed.  Add  data  cards  as 
necessary  to  contain  data. 


DATA  BLOCK  Q:  OMIT  IF  NPV  -  0  IN  BLOCK  B 

DESCRIPTION:  Variables  Co  be  changed  in  Che  paramecric  study 
FORMAT  AND  EXAMPLE 

1  2  3  4  5  6  7  8 


IPV  NIPV 


8  FORMAT 


2110 


VPI 

VP2 

VP3 

VP4 

NOTE:  READ  ONE  SET  OF  DATA  FOR  EACH  OF  THE  NPV  PARAMETRIC  VARIABLES. 
NOTE:  TWO  OR  MORE  DATA  CARDS  ARE  READ  FOR  EACH  SET  OF  DATA. 

FIELD  CONTENTS 

1  IPV:  Global  variable  number  associated  with  this  parametric 

variable. 

2  NIPV:  Number  of  separace  values  of  this  variable  to  be  read  on 

the  next  card. 


FIELD 


CONTENTS 


1-8  VPI:  I-th  value  of  the  parameter.  I  »  1  is  the  nominal  value 

about  which  the  parametric  study  is  performed. 


REMARKS 


1)  More  than  eight  values  of  the  parameter  are  allowed.  Add  data 
cards  as  required  to  contain  the  data. 


DATA  BLOCK  R 


DESCRIPTION:  COPES  data  "END“  card. 

FORMAT  AND  EXAMPLE 
1 

END 
END 

FIELD  CONTENTS 

1  END:  The  word  "END"  in  columns  1-3. 

REMARKS 

1)  This  card  MUST  appear  at  the  end  of  the  COPES  data. 

2)  This  ends  the  COPES  input  data. 

3)  Data  for  the  user-supplied  subroutine  ANALIZ  follows  this 


6.0  SAMPLE  COMPUTER  SOLUTIONS 


Sample  solutions  are  presented  here  for  the  variois  COPES  program 
options  defined  by  the  parameter  NCALC  In  DATA  BLOCK  B*  The  ANALIZ 
subroutine  given  in  Section  1.1  was  used  to  produce  the  results  given 
here.  Note  that  Subroutine  ANALIZ  does  not  read  data,  so  only  COPES 
data  is  required  here.  In  the  usual  case  where  data  is  read  as  input  to 
the  analysis  routine,  this  data  would  follow  directly  after  the  COPES 
"END"  card. 

The  output  from  a  COPES  program  execution  includes  a  title  page 
followed  by  a  copy  of  the  input  data.  Then  the  required  program 
executions  are  performed  and  the  final  output  information  is  printed. 
Figures  2-5  contain  output  for  all  options  of  COPES  as  follc«s: 


NCALC 

FIGURE 

PROBLEM  SOLVED 

1 

2 

Analysis  only* 

2 

3 

Optimization  by  ADS. 

3 

4 

Optimization  using  approximation 

techniques. 

4 

5 

Parametric  study. 

Note  that  data  for  the  COPES  program  is  read  depending  on  the  value 
of  the  parameters  NDV,  NAPRX  and  NPV  in  DATA  BLOCK  B.  The  actual 
program  execution  is  determined  by  the  value  of  NCALC.  Therefore,  data 
may  be  read  for  all  program  options  even  though  only  one  option  of  the 
program  is  to  be  executed.  Figure  6  is  a  copy  of  the  combined  data  for 
all  the  previous  examples.  This  data  may  be  used  for  any  program  option 
simply  by  changing  the  value  of  NCALC. 


n  n  n  n  o 


CCCCCCC  0000000 

0  0 
0  0 
0  0 
0  0 
0  0 
CCCCCCC  0000000 


ppppppp 

p  p 

p  p 

ppppppp 

p 

p 

p 


EEEEEEE 

E 

E 

EEEE 

E 

E 

EEEEEEE 


COPES/ADS 
CONTROL  PROGRAM 
FOR 

ENGINEERING  SYNTHE 
USING  THE 

ADS  OPTIMIZATION  PR 


TITLE 

THREE-BAR  TRUSS 

CARD  IMAGES  OF  INPUT  DATA 
CARD  IMAGE 

1)  THREE-BAR  TRUSS 

2)  1, 

3 )  END 


SSSSSSS 

S 

S 

SSSSSSS 

S 

S 

SSSSSSS 


SIS 


GRAM 


Figure  2:  COPES  Example,  Analysis  Only 


TITLE: 

THREE-BAR  TRUSS 
CONTROL  PARAMETERS; 

CALCULATION  CONTROL,  NCALC  -  1 

NUMBER  OF  GLOBAL  DESIGN  VARIABLES,  NDV  -  0 

INPUT  PARAMETER  FOR  APPROX.  OPT.,  NAPRX  -  0 

NUMBER  OF  PARAMETRIC  VARIABLES,  NPV  -  0 

INPUT  INFORMATION  PRINT  CODE,  IPNPUT  ■  0 

DEBUG  PRINT  CODE,  IPDBG  -  0 


CALCULATION  CONTROL,  NCALC 
VALUE  MEANING 

1  SINGLE  ANALYSIS 

2  OPTIMIZATION 

3  APPROXIMATE  OPTIMIZATION 

4  PARAMETRIC  STUDY 


*  *  ESTIMATED  DATA  STORAGE  REQUIREMENTS 


INPUT 

REAL 

EXECUTION 

AVAILABLE 

INPUT 

INTEGER 

EXECUTION 

7 

7 

5000 

1 

1 

THREE-BAR  TRUSS  OUTPUT 

Ai  «  0 . IOOOOE+0 l 
A2  -  0 . lOOOOE+Ql 
A3  -  0. IOOOOE+OI 

STRESSES 

SIGH  -  0.L4I42E+05 
SIG21  -  0.82843E+04 
SIG3I  — 0.58579E+04 
SIGI2  — 0.58579E+04 
SIG22  -  0.82843E+04 
SIG32  -  0.14142E+05 


VOLUME  *  0 . 38284E+02 
PROGRAM  CALLS  TO  ANALIZ 


ICALC  CALLS 

1  1 

2  I 

3  1 


AVAILABLE 

1000 


Figure  2  Concluded:  COPES  Example,  Analysis  Only 


CCCCCCC  0000000  PPPFPPP  EEEEEEE  SSSSSSS 

C  0  0  P  P  E  .  S 

C  0  0  P  P  E  S 

C  00  PPPPPPP  EEEE  SSSSSSS 

C  OOP  E  S 

C  OOP  E  S 

CCCCCCC  0000000  P  EEEEEEE  SSSSSSS 


COPES/ADS 
CONTROL  PROGRAM 
FOR 

ENGINEERING  SYNTHESIS 

USING  THE 

ADS  OPTIMIZATION  PROGRAM 


TITLE 

THREE-BAR  TRUSS 

CARD  IMAGES  OF  INPUT  DATA 

CARD  IMAGE 

1)  THREE-BAR  TRUSS 

2)  $  BLOCK  B 

3)  2,2 

4)  $  BLOCK  C 

5)  0,5, 7,0 

6)  $  BLOCK  D 

7)  3, 4, -I. 

8)  $  BLOCKS  E  AND  F  NOT  NEEDED 

9)  $  BLOCK  G 

10)  .01,1.0+20 

11)  .01,1.0+20 

12)  $  BLOCK  H 

13)  1,1 

14)  2,2 

15)  1,3 

16)  $  BLOCK  I 

17)  1 

18)  $  BLOCK  J 

19)  5,10 

20)  -15000. ,, 20000. 

21)  $  BLOCKS  K-0  NOT  NEEDED 

22)  $  BLOCK  R 

23)  END 


Figure  3:  COPES  Example,  Optimization 


TITLES 

THREE-BAR  TRUSS 


CONTROL  PARAMETERS; 

CALCULATION  CONTROL,  NCALC  -  2 
NUMBER  OF  GLOBAL  DESIGN  VARIABLES,  NDV  -  2 
INPUT  PARAMETER  FOR  APPROX.  OPT.,  NAPRX  -  0 
NUMBER  OF  PARAMETRIC  VARIABLES,  NPV  -  0 
INPUT  INFORMATION  PRINT  CODE,  IPNPUT  -  0 
DEBUG  PRINT  CODE,  IPDBG  -  0 


CALCULATION  CONTROL,  NCALC 
VALUE  MEANING 

1  SINGLE  ANALYSIS 

2  OPTIMIZATION 

3  APPROXIMATE  OPTIMIZATION 

4  PARAMETRIC  STUDY 

*  *  OPTIMIZATION  INFORMATION 
OPTIMIZATION  CONTROL  PARAMETERS 

ISTRAT  IOPT  IONED  IPRINT  NSCAL  IOVRDR  IOVRDI 
0  5  7  0  0  0  0 

GLOBAL  VARIABLE  NUMBER  OF  OBJECTIVE  -  4 

MULTIPLIER  (NEGATIVE  INDICATES  MINIMIZATION)  -  -0 .  lOOOE-i  "U 

DESIGN  VARIABLE  INFORMATION 

NON-ZERO  INITIAL  VALUE  WILL  OVER-RIDE  MODULE  INPUT 
D.  V.  LOWER  UPPER  INITIAL 

NO.  BOUND  BOUND  VALUE 

1  0 . 10000E-0 1  0. 11Q00E+I6  0.0 

2  O.IOOOQE-OL  0 . 1 1000E+I 6  0.0 

DESIGN  VARIABLES 

D.  V.  GLOBAL  MULTIPLYING 

ID  NO.  VAR.  NO.  FACTOR 

1  1  L  0 . 1 OOOOE+Q I 

22  2  0 . 1000QE+Q I 

31  3  0 . iOOQOE+O I 

CONSTRAINT  INFORMATION 

THERE  ARE  1  CONSTRAINT  SETS 

GLOBAL  GLOBAL  TYPE  LOWER  NORM.  UPPER  NORM. 

ID  VAR.  1  VAR.  2  ID  BOUND  FACTOR  BOUND  FACTOR 

L  5  10  0  -0 . 1 500E+0  5  0.1500E+05  0.2Q00E+05  0.2000E+05 

TOTAL  NUMBER  OF  CONSTRAINED  PARAMETERS  >  6 


Figure  3  Continued:  COPES  Example,  Optimization 


*  *  ESTIMATED  DATA  STORAGE  REQUIREMENTS 

REAL  INTEGER 

INPUT  EXECUTION  AVAILABLE  INPUT  EXECUTION  AVAILABLE 

42  192  5000  25  59  1000 


THREE-BAR  TRUSS  OUTPUT 

A1  -  0.10000E+01 
A2  -  0.10000E+01 
A3  -  0.10000E+01 

STRESSES 

SIGH  -  0. 14142E+05 
SIG21  -  0.82843E+04 
SIG31  — 0.58579E+04 
SIG12  — 0.58579E+04 
SIG22  -  0 . 82843E+04 
SIG32  -  0 . 14 142E+05 


VOLUME  -  0.38284E+02 


OPTIMIZATION  RESULTS 
OBJECTIVE  FUNCTION 

GLOBAL  LOCATION  4  FUNCTION  VALUE  0.26396E+02 

DESIGN  VARIABLES 


D.  V. 

GLOBAL 

LOWER 

UPPER 

ID 

NO. 

VAR.  NO. 

BOUND 

VALUE 

BOUND 

1 

1 

1 

0. 10000E-01 

0. 78474E+00 

0. 11000E+1 

2 

2 

2 

0. 10000E-01 

0 . 42005E+00 

0. 11000E  +  1 

3 

1 

3 

0. 10000E-01 

0. 78474E+00 

0. 11000E+1 

DESIGN 

CONSTRAINTS 

GLOBAL 

LOWER 

UPPER 

ID 

VAR.  NO. 

BOUND 

VALUE 

BOUND 

l  5  -0 . 15000E+0  5  0.199  9  6E+0  5  0.20000E+05 

3  6  -0 . 1 5000E+05  Q.14506E+05  0.20000E+05 

5  7  -0 . 15000E+05  -0.54903E+04  0.20000E+05 

7  8  0.1 5000E+05  -0.54903E+04  0.20000E+05 

9  9  -0 . 15000E+05  0.14506E+05  0.20000E+05 

11  10  -0 . 15000E+05  0 . 19996E+05  0.20000E+05 


Figure  3  Continued:  COPES  Example,  Optimization 


THREE-BAR  TRUSS  OUTPUT 


A1  -  0.7847 4E+00 
A2  -  0.42005E+00 
A3  -  0. 78474E+00 

STRESSES 

SIGH  -  0 . 19996E+05 
SIG21  -  0 . 14506E+05 
SIG3I  — 0.54903E+04 
SIG12  — 0.54903E+04 
SIG22  -  0 . 14506E+05 
SIG32  «  0 . 19996E+05 

VOLUME  -  0.26396E+02 


PROGRAM  CALLS  TO  ANALIZ 


ICALC  CALLS 

1  I 

2  21 

3  2 


Figure  3  Concluded:  COPES  Example, 


Optimization 


u  u  o  u  o 
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CCCCCCC 


0000000 
0  0 
0  0 
0  0 
0  0 
0  0 
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p  p 

p  p 

ppppppp 
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COPES/ADS 
CONTROL  PROGRAM 
FOR 

ENGINEERING  SYNTHE 
USING  THE 

ADS  OPTIMIZATION  PR 

TITLE 


THREE-BAR  TRUSS 


SSSSSSS 

S 

s 

sssssss 

s 

s 

sssssss 


SIS 


GRAM 


Figure  4:  COPES  Example,  Optimization  Using  Approximation  Techniques 


CARD  IMAGES  OF  INPUT  DATA 


CARD  IMAGE 

1)  THREE-BAR  TRUSS 

2)  $  BLOCK  B 

3)  3,2,1 

4)  $  BLOCK  C 

5)  0,5, 7,0 

6)  $  BLOCK  D 

7)  3, 4,-1. 

8)  $  BLOCKS  E  AND  F  NOT  NEEDED 

9)  $  BLOCK  G 
10)  .01,1.0+20 
11)  .01,1.0+20 

12)  $  BLOCK  H 

13)  1,1 

14)  2,2 

15)  1,3 

16)  $  BLOCK  I 

17)  l 

18)  $  BLOCK  J 

19)  5,10 

20)  -15000. , ,20000. 

21)  $  BLOCK  K 

22,  4,4 

23)  0 

24)  $  BLOCK  L 

25)  .2, .2 

26)  $  BLOCK  M 

27)  0. 

28)  $  BLOCK  N 

29)  .8, .6 

30)  .4,1. 

31)  1 .  ,  .  2 

32)  .5, .5 

33)  $  BLOCKS  0-S  NOT  NEEDED 

34)  $  BLOCK  R 

35)  END 


TITLE: 

THREE-BAR  TRUSS 


CONTROL  PARAMETERS; 

CALCULATION  CONTROL,  NCALC  -  3 
NUMBER  OF  GLOBAL  DESIGN  VARIABLES,  NDV  *  2 
INPUT  PARAMETER  FOR  APPROX.  OPT.,  NAPRX  -  1 
NUMBER  OF  PARAMETRIC  VARIABLES,  NPV  -  0 
INPUT  INFORMATION  PRINT  CODE,  IPNPUT  -  0 
DEBUG  PRINT  CODE,  IPDBG  -  0 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


■Ik  ,1m  A.  V. 


’**  '•*  %'  •  '  •  '  *,* 


CALCULATION  CONTROL,  NCALC 
VALUE  MEANING 

L  SINGLE  ANALYSIS 

2  OPTIMIZATION 

3  APPROXIMATE  OPTIMIZATION 

4  PARAMETRIC  STUDY 


*  *  OPTIMIZATION  INFORMATION 
OPTIMIZATION  CONTROL  PARAMETERS 

ISTRAT  IOPT  IONED  IPRINT  NSC AL  IOVRDR  IOVRDI 
0  5  7  0  0  0  0 

GLOBAL  VARIABLE  NUMBER  OF  OBJECTIVE  -  4 

MULTIPLIER  (NEGATIVE  INDICATES  MINIMIZATION)  -  -0.1000E+01 

DESIGN  VARIABLE  INFORMATION 

NON-ZERO  INITIAL  VALUE  WILL  OVER-RIDE  MODULE  INPUT 
D.  V.  LOWER  UPPER  INITIAL 

NO.  BOUND  BOUND  VALUE 

1  0.L0000E-01  0 . I IOOOE+I 6  0.0 

2  0. 10000E-01  0. 11000E+16  0.0 


DESIGN  VARIABLES 

D.  V.  GLOBAL  MULTIPLYING 

ID  NO.  VAR.  NO.  FACTOR 

II  I  0 . I0000E+0 1 

22  2  0 . 100Q0E+0 l 

31  3  0. 1Q000E+01 


CONSTRAINT  INFORMATION 
THERE  ARE  l  CONSTRAINT  SETS 

GLOBAL  GLOBAL  TYPE  LOWER  NORM.  UPPER  NORM. 

ID  VAR.  1  VAR.  2  ID  BOUND  FACTOR  BOUND  FACTOR 

1  5  10  0  -0 . 1 5 OOE+O 5  0.I500E+05  0.2000E+05  0.2000E+05 

TOTAL  NUMBER  OF  CONSTRAINED  PARAMETERS  -  6 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


*  *  APPROXIMATE  ANALYS IS /OPTIMIZATION  INFORMATION 


PRINT  CONTROL,  IPAPRX  -  1 
NUMBER  OF  INPUT  X-VECTORS,  NXI  -  4 
NUMBER  OF  INPUT  X-F  PAIRS,  NXFI  -  0 
X-VECTOR  FROM  ANALIZ ,  NXA  -  0 
NOMINAL  DESIGN,  INOM  -  0 
READ  UNIT  FOR  X-VECTORS,  ISCRX  -  5 
READ  UNIT  FOR  X-F  PAIRS,  ISCRXF  -  5 

MINIMUM  APPROXIMATING  CYCLES,  KMIN  -  1 
MAXIMUM  APPROXIMATING  CYCLES,  KMAX  -  12 
MAXIMUM  DESIGNS  USED  IN  FIT,  NPMAX  -  10 
NOMINAL  DESIGN  PARAMETER,  JNOM  -  10 
TAYLER  SERIES  I.D.  CODE,  MAXTRM  -  3 


DELTA-X  BOUNDS  FOR  APPROXIMATE  OPTIMIZATION 
0.2000E+00  0.2000E+00 

MULTIPLIER  ON  DELX,  XFACT1  -  0.1500E+01 

MULTIPLIER  ON  DELX,  XFACT2  -  0.2000E+01 


GLOBAL  LOCATIONS  OF  FUNCTIONS 

4  5  6  7  8  9  10 


X-VECTORS  INPUT  FROM  UNIT  5 

NUMBER  1  DESIGN  1 

0.8000E+00  0.6000E+00 


NUMBER  2  DESIGN  2 

0.4000E+00  0.1000E+01 


NUMBER  3  DESIGN  3 

0.1000E+01  0.2000E+00 


NUMBER  4  DESIGN  4 

0 . 500QE+00  0.5000E+00 


*  *  ESTIMATED  DATA  STORAGE  REQUIREMENTS 

REAL  INTEGER 

INPUT  EXECUTION  AVAILABLE  INPUT  EXECUTION  AVAILABLE 

42  255  5000  32  74  1000 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximati 

Technique  s 


APPROXIMATE  OPTIMIZATION  ITERATION  HISTORY 
APPROXIMATING  FUNCTION  I  IS  THE  OBJECTIVE 


BEGIN  ITERATION  NUMBER  I 

NOMINAL  DESIGN  NUMBER  «  I 
X-VECTOR 

0 . 80000E+Q0  0 . 60000E+00 

FUNCTION  VALUES 

0.28627E+02  0.18566E+05  0.12132E+05  -0.64340E+04  -0.64340E+04 

0.12132E+05  0 . 18566E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION,  COND  -  0.23630E-01 

RESULTS  OF  APPROXIMATE  OPTIMIZATION 

DELTA-X  VECTOR 
-0 . 522I2E-0I  -0 . 62160E-01 

X-VECTOR 

0.7477  9E+00  0.53784E+00 

APPROXIMATE  FUNCTION  VALUES 

0.26529E+02  0 . 20000E+0 5  0.13474E+05  -0.65257E+04  -0.65257E+04 

0 . 1 34  74E+05  0.20000E+05 

PRECISE  FUNCTION  VALUES 

0.26529E+02  0.20002E+05  0.13259E+05  -0.67433E+04  -0.67433E+04 

0 . 13259E+05  0.20002E+05 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


begin  iteration  number  2 

NOMINAL  DESIGN  NUMBER  -  5 


X-VECTOR 

0.7477  9E+00  0.53784E+00 


FUNCTION  VALUES 

0 • 26529E+02  0.20002E+05  0.13259E+05 

0. 13259E+05  0.20002E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION, 

RESULTS  OF  APPROXIMATE  OPTIMIZATION 


-0.67433E+04  -i 
COND  -  0.76112E 


DELTA-X  VECTOR 

0.64294E-01  -0.30000E+00 

X-VECTOR 

0.81208E+00  0 . 23  784E+00 

APPROXIMATE  FUNCTION  VALUES 
0.25348E+02  0.19996E+05 

0 • 1 5840E+Q5  0.19996E+05 

PRECISE  FUNCTION  VALUES 
0 . 25348E+02  0.2I02LE+05 

0.174 15E+05  0 . 2102IE+05 


0 . 15840E+05  -0.4I580E+04 

0.I7415E+05  -0.36066E+04 


.  6  7  43  3E+04 

03 


0.41 560E+04 

0 . 36066E+04 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


BEGIN  ITERATION  NUMBER  3 


NOMINAL  DESIGN  NUMBER  -  6 

X-VECTOR 

0 . 81 208E+00  0 . 23  784E+00 

FUNCTION  VALUES 

0.25348E+02  0.21021E+05  0.L7415E+05  -0.36066E+04  -0.36066E+04 

0.17  4 15E+05  0. 21021E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION,  COND  -  0.13815E-06 


RESULTS  OF  APPROXIMATE  OPTIMIZATION 

DELTA-X  VECTOR 
-0.22536E-01  0.15137E+00 

X-VECTOR 

0. 78955E+00  0.38921E+00 

APPROXIMATE  FUNCTION  VALUES 

0.26224E+02  0.20000E+05  0.14977E+05  -0.50230E+04  -0.50230E+04 
0.14977E+05  0.20000E+05 

PRECISE  FUNCTION  VALUES 

0.26224E+02  0.20128E+05  0.14926E+05  -0.52027E+04  -0.52027E+04 
0.14926E+05  0.20128E+05 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


BEGIN  ITERATION  NUMBER 


4 


NOMINAL  DESIGN  NUMBER  -  7 

X* V  EC  TOR 

0.78955E+00  0.38921E+00 

FUNCTION  VALUES 

0. 26224E+02  0.20128E+05  0 . 14 92 6E+05  -0.52027E+04  -0.52027E+04 

0 . 14926E+05  0.20128E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION,  COND  -  0.24905E-05 


RESULTS  OF  APPROXIMATE  OPTIMIZATION 


OPTIMIZATION  HAS  PRODUCED  AN  X-VECTOR  WHICH  IS  THE  SAME  AS  A  PREVIOUS 
DESIGN 

DELTA-X  VECTOR 

0.0  0.0 


X-VECTOR 

0.78955E+00  0.38921E+00 

THE  FOLLOWING  DESIGN  IS  NOT  THE  APPROXIMATE  OPTIMUM 

DELTA-X  VECTOR 

0 . 80000E-02  0.77921 E-Q2 

X-VECTOR 

0.797  55E+00  0. 39701E+00 

APPROXIMATE  FUNCTION  VALUES 

0 . 2652 8E+02  0.19871E+05  0.14711E+05  -0.5L594E+04  -0.51594E+04 

0.147 11E+05  0.1987 1E+05 

PRECISE  FUNCTION  VALUES 

0 . 26528E+02  0.19897E+05  0.14717E+05  -0.51801E+04  -0.51801E+04 

0.147 17E+05  0. 19897E+05 
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BEGIN  ITERATION  NUMBER 


5 


NOMINAL  DESIGN  NUMBER  -  7 

X-VECTOR 

0.78955E+00  0.38921E+00 

FUNCTION  VALUES 

0.26224E+02  0.20L28E+05  0.14926E+05  -0.52027E+04  -0.52027E+04 

0 . 14926E+05  0.20128E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION,  COND  -  0.25034E-05 


RESULTS  OF  APPROXIMATE  OPTIMIZATION 

DELTA-X  VECTOR 

0.47157E-02  0.16863E-02 

X-VECTOR 

0 . 7942  6E+00  0.39090E+00 

APPROXIMATE  FUNCTION  VALUES 

0.2637 4E+Q2  0.20000E+05  0.I4845E+05  -0.51550E+04  -0.51550E+04 
0 . I 4845E+05  0.20000E+05 

PRECISE  FUNCTION  VALUES 

0.2637 4E+02  0.200L4E+05  0.I4847E+05  -0.5I668E+04  -0.51668E+04 
0. I4847E+05  0.20014E+05 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


BEGIN  ITERATION  NUMBER 


6 


NOMINAL  DESIGN  NUMBER  -  9 

X-VECTOR 

0.79426E+00  0.39090E+00 

FUNCTION  VALUES 

0.26374E+02  0.20014E+05  0.14847E+05  -0.5L668E+04  -0.51668E+04 

0. 14847E+05  0.20014E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION,  COND  •  0.23234E-05 


RESULTS  OF  APPROXIMATE  OPTIMIZATION 


OPTIMIZATION  HAS  PRODUCED  AN  X-VECTOR  WHICH  IS  THE  SAME  AS  A  PREVIOUS 
DESIGN 


DELTA-X  VECTOR 

0.0  0.0 


X-VECTOR 

0 . 7942  6E+00  0.39090E+00 

THE  FOLLOWING  DESIGN  IS  NOT  THE  APPROXIMATE  OPTIMUM 

DELTA-X  VECTOR 

0 . 80000E-02  0.7  809  0E-02 

X-VECTOR 

0 . 80226E+00  0.3987IE+00 

APPROXIMATE  FUNCTION  VALUES 

0.2667 8E+02  0.L9758E+05  0.14634E+05  -0.5I241E+04  -0.51241E+04 
0. 14634E+05  0 . 197 58E+05 

PRECISE  FUNCTION  VALUES 

0.2667 8E+02  0.19785E+05  0.14640E+05  -0.51448E+04  -0.51448E+04 
0. 14640E+05  0.197  85E+05 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


BEGIN  ITERATION  NUMBER 


7 


NOMINAL  DESIGN  NUMBER  -  9 

X— VECTOR 

0.79426E+00  0.39090E+00 

FUNCTION  VALUES 

0.26374E+02  0.20014E+05  0.14847E+05  -0.51668E+04  -0.51668E+04 

0.14847E+05  0.20014E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION,  COND  ■  0.23355E-05 
SIDE  CONSTRAINTS  ON  APPROXIMATE  OPTIMIZATION 
LOWER  BOUNDS 

-O.IOOOOE+OO  -O.IOOOOE+OO 
UPPER  BOUNDS 

0 . lOOOOE+OO  O.IOOOOE+OO 


RESULTS  OF  APPROXIMATE  OPTIMIZATION 

DELTA-X  VECTOR 

0.0  0.0 


X-VECTOR 

0 . 79426E+00  0.39090E+00 

APPROXIMATE  FUNCTION  VALUES 

0.26374E+02  0.20014E+05  0.14847E+05  -0.51668E+04  -0.5I668E+04 

0 . I4847E+05  0.200L4E+05 


TWO  CONSECUTIVE  APPROXIMATE  OPTIMIZATIONS  HAVE  PRODUCED  THE  SAME  DESIGN 
OPTIMIZATION  TERMINATED 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 


PINAL  RESULT  OF  APPROXIMATE  OPTIMIZATION 

NOMINAL  DESIGN  NUMBER  -  9 

X-VECTOR 

0.79426E+00  0.39090E+00 

FUNCTION  VALUES 

0.26374E+02  0.20014E+05  0.14847E+05  -0.51668E+04  -0.51668E+04 

0.14847E+05  0.200I4E+05 

CONDITIONING  NUMBER  FOR  APPROXIMATION,  COND  -  0.23355E-05 


RESULTS  OF  APPROXIMATE  ANALYSIS /OPTIMIZATION 
TITLE 

THREE-BAR  TRUSS 

GLOBAL  LOCATIONS  OF  INDEPENDENT  X-VARIABLE S 
1  2 


GLOBAL  LOCATIONS  OF  FUNCTIONS,  F(X) 

4  5  6  7  8  9  10 

l 

APPROXIMATION  IS  BASED  ON  10  DESIGNS 

NOMINAL  DESIGN  IS  DESIGN  NUMBER  9 

VALUES  OF  X-V ARIABLES 

0.7943E+00  0 . 3909E+00 

VALUES  OF  FUNCTIONS,  F(X) 

0.2637E+02  0.2001E+05  0.1485E+05  -0.5167E+04  -0.5167E+04 

0 . 1485E+05  0.2001E+05 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 
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COEFFICIENTS  OF  TAYLOR  SERIES  EXPANSION 

PARAMETER  I  -  GLOBAL  VARIABLE  4 

LINEAR  TERMS,  DEL  F 

0. 2828E+02  0.1000E+02 

NON-LINEAR  TERMS,  H,  BEGINNING  WITH  DIAGONAL  ELEMENT 
ROW  1 

0. 1767E-02  0.9755E-03 

ROW  2 

0.3984E-03 


PARAMETER  2  -  GLOBAL  VARIABLE  5 

LINEAR  TERMS,  DEL  F 

-0.2412E+05  -0.8268E+04 

NON-LINEAR  TERMS,  H,  BEGINNING  WITH  DIAGONAL  ELEMENT 
ROW  I 

0. 4862E+05  0.4102E+04 

ROW  2 

0.I735E+05 


PARAMETER  3  -  GLOBAL  VARIABLE  6 

LINEAR  TERMS,  DEL  F 

-0.1129E+05  -0.1602E+05 

NON-LINEAR  TERMS,  H,  BEGINNING  WITH  DIAGONAL  ELEMENT 
ROW  1 

0.I633E+05  0.2327E+05 

ROW  2 

0.3I3IE+05 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniques 
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PARAMETER  4  -  GLOBAL  VARIABLE  7 

LINEAR  TERMS,  DEL  F 

0 . 1283E+05  -0. 7753E+04 

NON-LINEAR  TERMS,  H,  BEGINNING  WITH  DIAGONAL  ELEMENT 
ROW  L 

-0 . 3229E+05  Q.19I7E+05 

ROW  2 

0. 1396E+05 


PARAMETER  5  -  GLOBAL  VARIABLE  8 

LINEAR  TERMS,  DEL  F 

0. 1283E+05  -0. 7753E+04 

NON-LINEAR  TERMS,  H,  BEGINNING  WITH  DIAGONAL  ELEMENT 
ROW  L 

-0 . 3229E+05  0.19I7E+05 

ROW  2 

0. 1396E+05 


PARAMETER  6  -  GLOBAL  VARIABLE 

LINEAR  TERMS,  DEL  F 

-0. 1129E+05  -0 . I602E+05 


NON-LINEAR  TERMS,  H,  BEGINNING  WITH  DIAGONAL  ELEMENT 
ROW  I 

0.1633E+05  0.2327E+05 

ROW  2 

0.313 1E+05 


Figure  4  Continued:  COPES  Example,  Optimization  Using  Approximation 

Techniqu  es 
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PARAMETER 


7  -  GLOBAL  VARIABLE  10 


LINEAR  TERMS,  DEL  F 

-0 . 24 12E+05  -0.8268E+04 


NON-LINEAR  TERMS,  H,  BEGINNING  WITH  DIAGONAL  ELEMENT 
ROW  1 

0.4862E+05  0.4102E+04 


ROW  2 

0. 1735E+05 


OPTIMIZATION  RESULTS 


OBJECTIVE  FUNCTION 


GLOBAL 

LOCATION 

4 

FUNCTION  VALUE  0. 

26374E+02 

DESIGN 

VARIABLES 

D.  V. 

GLOBAL 

LOWER 

UPPER 

ID 

NO. 

VAR.  NO 

.  BOUND 

VALUE 

BOUND 

1 

1 

1 

-0 . 5 1448E+04 

0.79426E+00 

0. 11000E+1 

2 

2 

2 

0. 14640E+05 

0 . 39090E+00 

0. 11000E+1 

3 

1 

3 

-0. 51448E+04 

0. 79426E+00 

0. 11000E+1 

DESIGN  CONSTRAINTS 


ID 

I 

3 

5 

7 

9 

11 


GLOBAL 

LOWER 

VAR.  NO. 

BOUND 

5 

-0. 15000E+05 

6 

-0. 15000E+05 

7 

-0. 15000E+05 

8 

-0. 15000E+05 

9 

-0. 15000E+05 

10 

-0. 15000E+05 

VALUE 

0.20014E+05 
0. 14847E+05 
-0. 51668E+04 
-0. 51668E+04 
0. 14847E+05 
0. 20014E+05 


UPPER 

BOUND 

0. 20000E+05 
0. 20000E+05 
0. 20000E+05 
0. 20000E+05 
0. 20000E+05 
0. 20000E+05 


F igure 


4  Continued: 


COPES  Example,  Optimization  Using  Approximation 
Techniques 


THREE-BAR  TRUSS  OUTPUT 

A1  -  0. 79426E+Q0 
A2  -  0 • 39090E+00 
A3  -  0. 79426E+00 

STRESSES 

SIGH  -  0 . 200 14E+Q 3 
SIG21  -  0.14847E+05 
SIG3I  — 0.5I668E+04 
SIGI2  — 0.51668E+04 
SIG22  -  0.I4847E+05 
SIG32  -  0.20014E+05 

VOLUME  -  0. 26374E+02 


PROGRAM  CALLS  TO  ANALIZ 

ICALC  CALLS 
1  1 

2  10 

3  1 


Figure  4  Concluded:  COPES  Example, 

Techniques 


Optimization  Using  Approximati 
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OPTIMIZATION  PR 


TITLE 


THREE-BAR  TRUSS 


CARD  IMAGES  OF  INPUT  DATA 


CARD  IMAGE 

1)  THREE-BAR  TRUSS 

2)  $  BLOCK  B 

3)  4,, ,3 

4)  $  BLOCKS  C-0  NOT  NEEDED 

5)  $  BLOCK  P 

6)  7 

7)  4,5,6,7,8,9,10 

8)  $  BLOCK  Q 

9)  1,3 

10)  1., .8,1.2 

ID  2,3 

12)  l.,.8,l.2 

13)  3,5 

14)  1.,. 6, .8, 1.2, 1.4 

15)  $  BLOCK  R 

16)  END 


SSSSSSS 

S 

S 

SSSSSSS 

S 

S 

SSSSSSS 


SIS 


GRAM 


Figure  S:  COPES  Example,  Parametric  Study 


TITLE: 

THREE-BAR  TRUSS 
CONTROL  PARAMETERS; 

CALCULATION  CONTROL,  NCALC  -  4 

NUMBER  OF  GLOBAL  DESIGN  VARIABLES,  NDV  -  0 

INPUT  PARAMETER  FOR  APPROX.  OPT.,  NAPRX  -  0 

NUMBER  OF  PARAMETRIC  VARIABLES,  NPV  -  3 

INPUT  INFORMATION  PRINT  CODE,  IPNPUT  -  0 

DEBUG  PRINT  CODE,  IPDBG  -  0 

CALCULATION  CONTROL,  NCALC 
VALUE  MEANING 

1  SINGLE  ANALYSIS 

2  OPTIMIZATION 

3  APPROXIMATE  OPTIMIZATION 

4  PARAMETRIC  STUDY 

*  *  PARAMETRIC  INFORMATION 

PRINT  CONTROL,  IPRNTP  -  0 

NUMBER  OF  FUNCTIONS,  NPF  -  7 

GLOBAL  NUMBERS  ASSOCIATED  WITH  PARAMETRIC  FUNCTIONS 
4  5  6  7  8  9  10 

PARAMETERS 

GLOBAL  NOMINAL 

NUMBER  VARIABLE  VALUE  OFF-NOMINAL  VALUES 

1  I  0.10000E+01  0.8000E+00  0.1200E+01 

2  2  O.IOOOOE+OI  0 . 8000E+00  0.1200E+01 

3  3  O.IOOOOE+OI  0.6000E+00  0.8000E+00  0.I200E+01  0.1400E+01 


*  *  ESTIMATED  DATA  STORAGE  REQUIREMENTS 

REAL  INTEGER 

INPUT  EXECUTION  AVAILABLE  INPUT  EXECUTION  AVAILABLE 

18  25  5000  14  14  1000 


TITLE 

THREE-BAR  TRUSS 

NUMBER  OF  PARAMETRIC  VARIABLES,  NPV  -  3 

NUMBER  OF  FUNCTIONS,  NPF  -  7 

GLOBAL  LOCATIONS  OF  PARAMETRIC  VARIABLES 
12  3 

GLOBAL  LOCATIONS  OF  FUNCTIONS 

4  5  6  7  8  9  10 


Figure  5  Continued:  COPES  Example,  Parametric  Study 


NOMINAL  DESIGN  INFORMATION 

VALUES  OF  THE  PARAMETERS 

O.IOOOOE+Ol  0 . 10000E+01  O.IOOOOE+Ol 


VALUES  OF  THE  FUNCTIONS 

0 . 38284E+02  0.14142E+05  0.82843E+04  -0.58579E+04  -0.58579E+04 

0 . 82843E+04  0 . 14 I42E+05 

PARAMETRIC  STUDY  RESULTS 

GLOBAL  VARIABLE  1 


X 

F(X) 

0. 8000E+00 

0.3546E+02 
-0 . 6823E+04 

0. 1647E+05 

0. 77I9E+04 

0 . 9649E+04 
0. 1454E+05 

-0 . 6823E+04 

0. 1200E+01 

0.41 I1E+02 
-0.5132E+04 

0. 1239E+05 
0.8709E+04 

0. 7258E+04 
0. 1384E+05 

-0. 5I32E+04 

GLOBAL  VARIABLE 

2 

X 

F(X) 

0. 8000E+00 

0. 3628E+02 
-0 . 5308E+04 

0. 1469E+05 

0. 9384E+04 

0 . 9384E+04 
0. I469E+05 

-0. 5308E+04 

0. 1200E+01 

0 . 4028E+02 
-0 . 6292E+04 

0. 1371E+05 
0.7415E+04 

0.7415E+04 
0.  I371E+05 

-0.6292E+04 

GLOBAL  VARIABLE 

3 

X 

F  (X) 

0 . 60Q0E+00 

0 . 3  263E+02 
-0 . 8I68E+04 

0. L5I0E+Q5 

0. 1155E+05 

0 . 6931E+04 
0. 1972E+05 

-0.8I68E+04 

0 . 8000E+00 

0. 3546E+02 
-0.6823E+04 

0. 1454E+05 
0.9649E+04 

0. 7719E+04 
0. I647E+05 

-0.6823E+04 

0. 1200E+0I 

0.4II1E+02 
-0 . 5132E+04 

0.  I384E+05 

0. 7258E+04 

0.8709E+04 
0. 1239E+05 

-0. 5132E+04 

0. 1400E+01 

0.4394E+02 
-0. 4566E+Q4 

0. I361E+05 

0 . 6458E+04 

0. 9041E+04 
0. II02E+05 

-0 . 4566E+04 

PROGRAM  CALLS  TO  ANALIZ 


ICALC  CALLS 

1  I 

2  9 

3  0 


Figure  5  Concluded:  COPES  Exaaple,  Paraaetric  Study 


THREE-BAR  TRUSS 
$  BLOCK  B 

3. 2. 1.3 

$  BLOCK  C 
0,5, 7,0 
$  BLOCK  0 
3, 4,-1. 

$  BLOCKS  E  AND  F  NOT  NEEDED 
$  BLOCK  G 
.01,1.0+20 
.01,1.0+20 
$  BLOCK  H 
1,1 
2,2 

1.3 

$  BLOCK  I 
1 

$  BLOCK  J 
5,10 

-15000. , ,20000. 

$  BLOCK  K 

1.4  • 

0 

$  BLOCK  L 
.2,  .2 
$  BLOCK  M 
0. 

$  BLOCK  N 
.  8 ,  .  6 
•  4,1. 

1  .  ,  .2 

.5,  .  5 

$  BLOCK  0,  NO  DATA  BECAUSE  NXFI-0 
$  BLOCK  P 


4,5,6,7,8,9,10 
$  BLOCK  Q 

1.3 

1.  ,  .8,1.2 

2.3 

1., .8,1.2 

3,5 

1.,. 6, .8, 1.2, 1.4 

$  BLOCK  R 

END 


Figure  6:  COPES  Example,  Combined  Data  For  All  Options 
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APPENDIX  A 


ADS  PROGRAM  OPTIONS 


The  options  available  in  the  ADS  program  are  identified  here*  This 
information  is  taken  from  the  ADS  User's  manual  (ref.  4). 

A. 1  INTRODUCTION 

ADS  is  a  general  purpose  iumerical  optimization  program  containing 
a  wide  variety  of  algorithms.  The  problem  solved  is: 


Minimize  F(X) 

Subject  to; 

gj(X)  <  0 

j-i,m 

hk(X)  -  o 

k-1,1 

v!  a 

Xi  <  Xi  <  Xi 

i*l,n 

The  solution  of  this  general  problem  is  separated  into  three  basic 
levels : 

1.  STRATEGY  -  For  example,  Sequential  Unconstrained  Minimization  or 
Sequential  Linear  Programming. 

2.  OPTIMIZER  -  For  example.  Variable  Metric  methods  for  unconstrained 
minimization  or  the  Method  of  Feasible  Directions  for  constrained 
minimization. 

3.  ONE -DIMENSIONAL  SEARCH  -  For  example,  Golden  Section  or  Polynomial 
Interpolation. 

Additionally,  we  may  consider  another  component  to  be  problem 
formulation.  It  is  assumed  that  the  engineer  makes  every  effort  to 
formulate  the  problem  in  a  form  amenable  to  efficient  solution  by 
numerical  optimization.  This  aspect  is  perhaps  the  most  important 
ingredient  to  the  efficient  use  of  the  ADS  program  for  solution  of 
problems  of  practical  significance. 

By  choosing  the  Strategy,  Optimizer  and  One-Dimensional  Searcn,  the 
user  is  given  considerable  flexibility  in  creating  an  optimization 
program  which  works  well  for  a  given  class  of  design  problems. 

The  purpose  here  is  to  describe  the 

available  program  options.  Section sA.2 -A4  identify  the  available 
optimization  strategies,  optimizers  and  one-dimensional  search 
algorithms.  Section  A.5  gives  allowed  combinations  of  algorithms  and  a 
Quick-Reference  summary.  Section  A. 6  gives  a  brief  set  of  guidelines 
for  choosing  an  overall  optimization  algorithm. 


A.2  Strategy 

Table  A.1  lists  the  strategies  available.  The  parameter  ISTRAT 
will  be  seat  to  the  ADS  program  to  Identify  the  strategy  the  user  wants. 
The  1STRAT*0  option  would  indicate  that  control  should  transfer  directly 
to  the  optimizer.  This  would  be  the  case,  for  example,  when  using  the 
Method  of  Feasible  Directions  to  solve  constrained  optimization  problems 
because  the  optimizer  works  directly  with  the  constrained  problem.  On 
the  other  hand,  if  the  constrained  optimization  problem  is  to  be  solved 
by  creating  a  sequence  of  unconstrained  minimizations,  with  penalty 
functions  to  deal  with  constraints,  one  of  the  appropriate  strategies 
would  be  used. 


TABLE  A.l:  STRATEGY  OPTIONS 

ISTRAT  STRATEGY  TO  BE  USED 


0 

1 

2 

3 

4 

5 

6 

7 

8 


None.  Go  directly  to  the  optimizer. 

Sequential  unconstrained  minimization  using  the  exterior 
penalty  function  method. 

Sequential  unconstrained  minimization  using  the  linear 
extended  interior  penalty  function  method. 

Sequential  unconstrained  minimization  using  the  quadratic 
extended  interior  penalty  function  method. 

Sequential  unconstrained  minimization  using  the  cubic 
extended  interior  penalty  function  method. 

Augmented  Lagrange  Multiplier  method. 

Sequential  Linear  Programming. 

Method  of  Centers  (method  of  inscribed  hyperspheres) 
Sequential  Quadratic  Programming. 


A. 3  Optimizer 

Table  A.2  lists  the  optimizers  available.  10PT  is  the  parameter 
used  to  Indicate  the  optimizer  desired. 


TABLE  A.2:  OPTIMIZER  OPTIONS 

IOPT  OPTIMIZER  TO  BE  USED 

0  None.  Go  directly  to  the  one-dimensional  search.  This 
option  should  be  used  only  for  program  development. 

1  Fletcher-Reeves  algorithm  for  unconstrained  minimization 

2  Davidon-Fletcher-Powell  (DFP)  variable  metric  method  for 
unconstrained  minimization. 

3  Broydon-Fletcher-Goldfarb-Shanno  (BFGS)  variable  metric 
method  for  unconstrained  minimization. 

4  Method  of  Feasible  Directions  (MFD)  for  constrained 
minimization. 

5  Modified  Method  of  Feasible  Directions  for  constrained 
minimization. 
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la  choosing  the  optimizer  (as  well  as  strategy  and  one-dimensional 
search)  it  is  assumed  that  the  user  isicnow  ledge  able  enough  to  choose 
an  algorithm  consistent  with  the  problem  at  hand*  For  example,  a 
variable  metric  optimizer  would  not  be  used  to  solve  constrained 
problems  unless  a  strategy  is  used  to  create  the  equivalent 
unconstrained  minimization  task  via  some  form  of  penalty  function. 

A. 4  One-D lmens ional  Search 

Table  A.3  lists  the  one-dimensional  search  options  available  for 
unconstrained  and  constrained  problems.  Here  IONED  identifies  the 
algorithm  to  be  used. 

TABLE  A.3:  ONE-DIMENSIONAL  SEARCH  OPTIONS 
IONED  ONE-DIMENSIONAL  SEARCH  OPTION 

1  Find  the  minimum  of  an  unconstrained  function  using  the 
Golden  Section  method. 

2  Find  the  minimum  of  an  unconstrained  function  using  the 
Golden  Section  method  followed  by  polynomial  interpolation. 

3  Find  the  minimum  of  an  unconstrained  function  by  first 
finding  bounds  and  then  using  polynomial  interpolation. 

4  Find  the  minimum  of  an  unconstrained  function  by  polynomial 
interpolation/extrapolation  without  first  finding  bounds  on 
the  solution. 

3  Find  the  minimum  of  an  ‘constrained  function  using  the 
Golden  Section  method. 

6  Find  the  minimum  of  an  constrained  function  using  the 
Golden  Section  method  followed  by  polynomial  interpolation. 

7  Find  the  minimum  of  an  constrained  function  by  first 
finding  bounds  and  theniusing  polynomial  interpolation. 

8  Find  the  minimum  of  an  constrained  function  by  polynomial 
interpolation/extrapolation  without  first  finding  bounds  on 
the  solution. 


A. 5  Allowable  Combinations  of  Algorithms 

Not  all  combinations  of  strategy,  optimizer  and  one-dimensional 
search  are  meaningful.  For  example,  constrained  one-dimensional  search 
is  not  meaningful  when  minimizing  unconstrained  functions. 

Table  A.4  Identifies  the  combinations  of  algorithms  which  are 
available  in  the  ADS  program.  In  this  table,  an  X  is  used  to  denote  an 
acceptable  combination  of  strategy,  optimizer  and  one-dimensional 
search.  An  example  is  shown  by  the  heavy  line  on  the  table  which 
indicates  that  constrained  optimization  is  to  be  performed  by  the 
Augmented  Lagrange  Multiplier  Method  (ISTRAT*5),  using  the  BFGS 
optimizer  (I0PT»3)  and  polynomial  interpolation  with  bounds  for  the  one- 
dimensional  search  (I0NED*3).  From  the  table,  it  is  clear  that  a  lazge 
number  of  possible  combinations  of  algorithms  are  available. 


TABLE  A. 5:  QUICK  REFERENCE  TO  ADS  OPTIONS 


IOPT  OPTIMIZER 


STRATEGY  ISTRAT  IOPT  12345 


None 

0 

X 

X 

X 

X 

X 

SUMT,  Exterior 

1 

X 

X 

X 

0 

0 

SUMT,  Linear  Extended  Interior 

2 

X 

X 

X 

0 

0 

SUMT,  Quadratic  Extended  Interior 

3 

X 

X 

X 

0 

0 

SUMT,  Cubic  Extended  Interior 

4 

X 

X 

X 

0 

0 

Augmented  Lagrange  Multiplier  Meth 

.  5 

X 

X 

X 

0 

0 

Sequential  Linear  Programming 

6 

0 

0 

0 

X 

X 

Method  of  Centers 

7 

0 

0 

0 

X 

X 

Sequential  Quadratic  Programming 

8 

0 

0 

0 

X 

X 

ONE— DIMENSIOIIAL  SEARCH  . 

IONED 

Golden  Section  Method  l  X  X  X  0  0 
(olden  Section  +  Polynomial  2  X  X  X  0  0 
Polynomial  Interpolation  (bounded)  3  X  X  X  0  0 
Polynomial  Extrapolation  4  X  X  X  0*0 
Golden  Section  Method  5  0  0  0  X  X 
Golden  Section  +  Polynomial  6  0  0  0  X  X 
Polynomial  Interpolation  (bounded)  7  0  0  0  X  X 
Polynomial  Extrapolation  8  0  0  0  X  X 


NOTE:  An  X  denotes  an  allowed  combination  of  algorithms' 


A. 6  Choosing  An  Algorithm 

One  difficulty  with  a  program  such  as  ADS,  which  provides  numerous 
options,  is  that  of  picking  the  best  combination  of  algorithms  to  solve 
a  given  problem.  While  it  is  not  possible  to  provide  a  concise  set  of 
rules,  some  general  guidelines  are  offered  here  based  on  the  author's 
experience.  The  user  is  strongly  encouraged  to  try  many  different 
options  in  order  to  gain  familiarity  with  ADS  and  to  Improve  the 
probability  that  the  best  combination  of  algorithms  is  found  for  the 
particular  class  of  problems  being  solved. 


UNCONSTRAINED  FUNCTIONS  (NCON-O,  Side  Constraints  OK) 

ISTRAT-0 

Is  computer  storage  very  limited? 

Tes  -  IOPT-1.  Are  function  evaluations  expensive? 

Yes  -  Is  the  objective  known  to  be  approximately  quadratic? 

Yes  -  IONED-4 
No  -  IONED-3 
No  -  IONED-1  or  2 
No  -  Is  the  analysis  iterative? 

Yes  -  IQPT-3.  Are  function  evaluations  expensive? 

Yes  -  Is  the  objective  known  to  be  approximately  quadratic?  . 

Yes  -  IONED-4 
No  -  IONED-3 
No  -  IONED-1  or  2 

No  -  IOPT-2  or  3.  Are  function  evaluations  expensive? 

Yes  -  Is  the  objective  known  to  be  approximately  quadratic? 

Yes  -  IONED-4 
No  -  IONED-3 
No  -  IONED-1  or  2 


CONSTRAINED  FUNCTIONS 

Are  relative  minima  known  to  exist? 

Yes  -  ISTRAT-1,  IOPT-3.  Are  function  evaluations  expensive? 
Yes  -  IONED-3 
No  -  IONED-1  or  2 

No  -  Are  the  objective  and/or  constraints  highly  nonlinear? 
Yes  -  Are  function  evaluations  expensive? 

Yes  ~  ISTRAT-O,  IOPT-4,  IONED-7 
No  -  ISTRAT-2,  3  or  5,  IOPT-2  or  3,  IONED-1  or  2 
No  -  Is  the  design  expected  to  be  fully-constrained? 
(i.e.  NDV  active  constraints  at  the  optimum) 

Yes  -  ISTRAT-6 ,  IOPT-5,  IONED-6 
No  -  Is  the  analysis  iterative? 

Yes  ~  ISTRAT-O,  IOPT-4,  IONED-7  or 
ISTRAT-8,  IOPT-5,  IONED-7 
No  -  ISTRAT-O,  IOPT-5,  IONED-7  or 
ISTRAT-8,  IOPT-5,  IONED-7 


GENERAL  APPLICATIONS 


k. 


| 

l\* 


Often  little  is  known  about  the  nature  of  the  problem  being  solved. 
Based  on  experience  with  a  wide  variety  of  problems,  a  very  direct 
approach  is  given  here  for  using  ADS.  The  following  table  of  parameters 
is  offered  as  a  sequence  of  algorithms.  When  using  ADS  the  first  few 
times,  the  user  may  prefer  to  run  the  cases  given  here,  rather  than 
using  the  decision  approach  given  above.  It  is  assumed  here  that 
a  constrained  optimization  problem  is  being  solved.  If  the  problem  is 
unconstrained,  ISTRAT-0,  I0PT-3  and  IONED-2  or  3  is  recommended. 


ISTRAT 

IOPT 

10NED 

IPRINT 

8 

5 

7 

2200 

0 

5 

7 

2020 

0 

4 

7 

2020 

6 

5 

6 

2200 

5 

3 

3 

2200 

2 

3 

3 

2200 

1 

3 

3 

2200 

69 


*  **•  "  *  *  •  *  •  •  *  •  •  »  *  •  *  •  *  »  m  m  *  >  "  •  *  m  %  *»  *  ■  ®  a  *  «  "  >.  “  m  *•  m  - 


APPENDIX  B 


OVER-RIDING  ADS  DEFAULT  PARAMETERS 

Various  internal  parameters  are  defined  on  the  first  call  to  ADS 
which  work  well  for  the  "average"  optimization  task.  However,  it  is 
often  desirable  to  change  these  in  order  to  gain  maximum  utility  of  the 
program.  Various  real  and  integer  scalar  parameters  are  stored  in  arrays 
WK  and  IHX  respectively.  Those  which  the  user  may  wish  to  change  are 
listed  in  Tables  B.l  through  B.4,  together  with  their  default  values  and 
definitions.  If  the  user  wishes  to  change  any  of  these,  the  appropriate 
component  of  HK  or  IWK  is  simply  re-defined.  For  example,  if  the 
relative  convergence  criterion,  DELOB J,  is  to  be  changed  to  0.002,  this 
is  done  by  setting  location  12  in  WK  to  0.002,  because  WK(12)  contains 
the  value  of  DEL0BJ. 


TABLE  B.l:  URAL  PARAMETERS  STORED  IN  ARRAY  UK 

MODULES  WHERE  USED 


PARAMETER 

LOCATION  DEFAULT 

ISTRAT 

IOPT 

IONED 

ALAMDZ 

1 

0.0 

5 

- 

- 

BETAMC 

2 

0.0 

7 

- 

- 

CT1 

3 

-0.03 

- 

4,5 

- 

CTL 

A 

-0.005 

- 

4,5 

- 

CTLMIN 

5 

0.001 

- 

4,5 

- 

CTMIN 

6 

0.01 

- 

4,5 

- 

DAB  ALP2 

7 

0.0001 

- 

ALL 

— 

DABOBJ 

8 

ABS(F0)/10000 

ALL 

- 

- 

DABOBM 

9 

ABS(F0)/1000 

ALL 

— 

— 

DABSTR 

10 

ABS(F0)/10000 

ALL 

- 

DELALP3 

11 

0.005 

- 

- 

1.2, 5, 6 

DELOB J 

12 

0.001 

- 

ALL 

• 

DELOBM 

13 

0.01 

ALL 

— 

— 

delstr 

14 

0.001 

ALL 

- 

- 

DLOBJ1 

15 

0.1 

— 

ALL 

• 

DL0BJ2 

16 

1000.0 

- 

ALL 

- 

DX1 

17 

0.1 

- 

ALL 

- 

DX2 

18 

0.2 

- 

ALL 

- 

EPS PEN 

19 

-0.05 

2,3,4 

- 

EXTRAP 

20 

5.0 

- 

- 

ALL 

FDCH 

21 

0.01 

- 

ALL 

— 

FDCHM 

22 

0.001 

- 

ALL 

- 

GMULTZ 

23 

10.0 

8 

- 

- 

PSAIZ 

24 

0.95 

8 

- 

- 

RMULT 

25 

5.0 

1,5 

- 

— 

RMVLMZ 

26 

0.2 

6,7,8 

- 

- 

RP 

27 

10.0 

1,5 

- 

- 

RPMAX 

28 

1.0E+10 

1,5 

- 

- 

RPMULT 

29 

0.2 

1,5 

- 

— 

RPPMIN 

30 

1.0E-10 

2,3,4 

— 

- 

RPPRIM 

31 

100.0 

2,3,4 

- 

- 

SCFO 

32 

1.0 

ALL 

ALL 

ALL 

SCLMIN 

33 

0.001 

ALL 

ALL 

ALL 

STOL 

34 

0.001 

- 

4,5 

- 

THETAZ 

35 

0.1 

- 

4,5 

- 

XMULT 

36 

2.618034 

- 

- 

1, 2,3,5, 

ZRO 

37 

0.00001 

ALL 

ALL 

ALL 

1  If  IOPT-4,  CT— 0.1 

2  If  IONED-3  or  8,  DABALP-0.001 

3  If  IONED-3  or  8,  DELALP-0.05 


NOTE:  FO  is  the  magnituge  of  the  initial  objective  function. 


TABLE  B.2:  DEFINITIONS  OF  REAL  PARAMETERS  CONTAINED  IN  ARRAY  WK 
PARAMETER _ DEFINITION _ 

ALAMDZ  Initial  estimate  of  the  Lagrange  Multipliers  in  the  Augmented 
Lagrange  Multiplier  Method* 

BETAMC  Additional  steepest  descent  fraction  in  the  method  of  centers. 
After  moving  to  the  center  of  the  hypersphere,  a  steepest 
descent  move  is  made  equal  to  BETAMC  times  the  radius  of  the 
hyper sphere* 

CT  Constraint  tolerance  in  the  Method  of  Feasible  Directions  or 

the  Modified  Method  of  Feasible  Directions.  A  constraint  is 
active  If  its  aimer ical  value  is  more  positive  than  CT. 

CTL  Same  as  CT,  but  for  linear  constraints. 

CTLMIN  Same  as  CTMIN,  but  for  linear  constraints. 

CTMIN  Minimum  constraint  tolerance  for  nonlinear  constraints.  If  a 
constraint  is  more  positive  than  CTMIN,  it  is  considered  to 
be  violated. 

DABALP  Absolute  convergence  criteria  for  the  one-dimensional  search 
when  using  the  Golden  Section  method. 

DABOBJ  Maximum  absolute  change  in  the  objective  between  two 
consecutive  iterations  to  indicate  convergence  in 
optimization. 

DABOBM  Absolute  convergence  criterion  for  the  optimization  sub¬ 
problem  when  using  sequential  minimization  techniques. 

DABSTR  Same  as  DABOBJ,  but  used  at  the  strategy  level. 

DELALP  Relative  convergence  criteria  for  the  one-dimensional  search 
when  using  the  Golden  Section  method. 

DELOBJ  Maximum  relative  change  in  the  objective  between  two 
consecutive  iterations  to  indicate  convergence  in 
optimization. 

DELOBM  Relative  convergence  criterion  for  the  optimization  sub¬ 
problem  when  using  sequential  minimization  techniques. 

DELSTR  Same  as  DELOBJ,  but  used  at  the  strategy  level. 

DL0BJ1  Relative  change  in  the  objective  function  attempted  on  the 

first  optimization  iteration.  Used  to  estimate  initial  move 
in  the  one-dimensional  search.  Updated  as  the 
optimization  progresses. 

DL0BJ2  Absolute  change  in  the  objective  function  attempted  on  the 

first  optimization  iteration.  Used  to  estimate  initial  move 
in  the  one-dimensional  search.  Updated  as  the 
optimization  progresses. 

DX1  Maximum  relative  change  in  a  design  variable  attempted  on  the 
first  optimization  iteration.  Used  to  estimate  the  initial 
move  in  the  one-dimensional  search.  Updated  as  the 
optimization  progresses. 

DX2  Maximum  absolute  change  in  a  design  variable  attempted  on  the 
first  optimization  iteration.  Used  to  estimate  the  initial 
move  in  the  one-dimensional  search.  Updated  as  the 
optimization  progresses. 

EPSPEN  Initial  transition  point  for  extended  penalty  function 
methods.  Updated  as  the  optimization  progresses. 

EXTRAP  Maximum  multiplier  on  the  one-dimensional  search  parameter, 
ALPHA  in  the  one-dimensional  search  using  polynomial 
interpolation/extrapolation. 


TABLE  B.2  CONCLUDED:  DEFINITIONS  OF  REAL  PARAMETERS  CONTAINED  IN  ARRAY  HK 
PARAMETER  DEFINITION 

FDCH  Relative  finite  difference  step  when  calculating  gradients. 

FDCm  Minimum  absolute  value  of  the  finite  difference  step  when 

calculating  gradients.  This  prevents  too  small  a  step  when 
X(I)  is  near  zero. 

GMUL1Z  Initial  penalty  parameter  in  Sequential  Quadratic  programming. 

PSA2Z  Move  fraction  to  avoid  constraint  violations  in  Sequential 
Quadratic  Programming. 

KMULT  Penalty  function  nultiplier  for  the  exterior  penalty  function 
method.  Must  be  greater  than  1.0. 

RMVLMZ  Initial  relative  move  limit.  Used  to  set  the  move  limits  in 
sequential  linear  programming,  method  of  inscribed 
hyperspheres  and  sequential  quadratic  programming  as  a 
fraction  of  the  value  of  X(I),  I»1,NDV. 

RP  Initial  penalty  parameter  for  the  exterior  penalty  function 

method  or  the  Augmented  Lagrange  Multiplier  method. 

RPMAX  Maximum  value  of  RP  for  the  exterior  penalty  function  method 
or  the  Augmented  Lagrange  Multiplier  method. 

RPMULT  Multiplier  on  RP  for  consequtive  iterations. 

RRPMIN  Minimum  value  of  RPPRIM  to  indicate  convergence. 

RPPRIM  Initial  penalty  parameter  for  extended  interior  penalty 
function  methods. 

SCFO  The  user-supplied  value  of  the  scale  factor  for  the  objective 
function’if  the  default  or  calculated  value  is  to  be  over¬ 
ridden. 

.  SCLMIN  Minimum  numerical  value  of  any  scale  factor  allowed. 

STOL  Tolerance  on  the  components  of  the  calculated  search  direction 
to  Indicate  that  the  Kuhn-Tucker  conditions  are  satisfied. 

THETAZ  Nominal  value  of  the  push-off  factor  in  the  Method  of  Feasible 
Directions. 

XMULT  Multiplier  on  the  move  parameter,  ALPHA,  in  the  one¬ 
dimensional  search  to  find  bounds  on  the  solution. 

ZRO  Numerical  estimate  of  zero  on  the  computer.  Usually  the 

default  value  is  adequate.  If  a  computer  with  a  short  word 
length  is  used,  ZRO*1.0E-4  may  be  preferred. 


TABLE  B.3: 

INTEGER  PARAMETERS 

STORED  IN 

ARRAY 

IWK 

MODULES 

WHERE 

USED 

PARAMETER 

LOCATION 

DEFAULT 

ISTRAT 

I0PT 

IONED 

ICNDIR 

1 

NDV+1 

- 

ALL 

- 

ISCAL 

2 

1 

ALL 

ALL 

ALL 

UMAX 

3 

40 

- 

ALL 

- 

ITRMOP 

4 

3 

- 

1,2,3 

- 

ITRMST 

5 

2 

ALL 

- 

JONED 

6 

IONED 

8 

- 

- 

JTMAX 

7 

20 

ALL 

- 

- 

TABLE  B.4:  DEFINITIONS  OF  INTEGER  PARAMETERS  CONTAINED  IN  ARRAY  INK 
PARAMETER _ DEFINITION _ 

ICNDIR  Restart  parameter  for  conjugate  direction  and  variable  metric 
methods*  Unconstrained  minimization  is  restarted  with  a 
steepest  descent  direction  every  ICNDIR  iterations. 

ISCAL  Scaling  parameter.  If  ISCAL~0,  no  scaling  is  done.  If 

ISCAL-1,  the  design  variables,  objective  and  constraints  are 
scaled  automatically. 

UMAX  Maxioum  number  of  iterations  allowed  at  the  optimizer  level. 
ITRMOP  The  number  of  consecutive  iterations  for  which  the  absolute  or 
relative  convergence  criteria  must  be  met  to  indicate 
convergence  at  the  optimizer  level. 

ITRMST  The  number  of  consecutive  iterations  for  which  the  absolute  or 
relative  convergence  criteria  must  be  met  to  indicate 
convergence  at  the  strategy  level. 

JONED  The  one-dimensional  search  parameter  (IONED)  to  be  used  in  the 
Sequential  Quadratic  Programming  method  at  the  strategy  level. 
JTMAX  Maximum  number  of  iterations  allowed  at  the  strategy  level. 
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